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Abstract 

This  report  continues  the  investigation  of  Research  Report  No,  EM-/17 
on  the  applications  of  recently  developed  mathematical  techniques.  The  aim 
here  is  to  present  the  theory  of  Green's  f-ionctions  for  ordinary  differentisil 
equations  and  for  systems  of  ordinary  differential  equations.  This  theory 
can  "be  used  to  ohtain  the  spectral  representation  of  ordinary  differential 
eqimtions  and  to  solve  partial  differential  equations. 
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Introduction 


In  this  report  the  methods  previously  developed  will  be  used  to  discuss 


the  solutions  of  ordinary  linear  differential  equations.  The  main  problem  will 


be  to  solve  y 


(3.1)        Lu  =  f, 

where  L  is  a  differential  operator,  f  a  known  fiinction,  and  u  the  unknown  f\inction. 
The  general  theory  of  linear  equations  tells  us  that  there  are  two  important 
approaches  to  solving  (3.1).  One  approach  is  to  find  the  inverse  operator  to  L, 
that  is,  to  find  an  operator  L~  such  that  L"  L  is  the  identity  operator.  The 
other  approach  is  to  find  the  spectral  representation  of  L  by  studying  the  solu- 
tions of  the  equations 

Lu  =  \u, 

where  X  is  an  arbitrary  constant.  Of  course,  L  is  not  a  finite-dimensional 
bounded  operator,  so  the  theorems  of  matrix  theory  are  not  imiriediately  valid, 
We  shall  use  them,  however,  to  suggest  the  theorems  which  we  shall  try  to  prove. 

Both  approaches,  that  of  finding  the  inverse  operator  and  that  of  finding 
its  spectral  representation  will  be  used  in  discussing  differential  operators. 
Most  of  the  discussion  will  be  devoted  to  the  consideration  of  gjecial  examples 
in  order  to  reveal  the  techniques  and  the  possible  insults.  After  the  treatment 
of  each  group  of  special  examples,  there  will  be  a  section  giving  a  general  dis- 
cussion of  the  methods  employed  in  that  particular  type  of  problem. 

The  Linear  Vector  Space  and  its  Operators 

The  linear  vector  space  that  will  be  used  is  the  function  space  of  all  real 
functions  which  are  Lebesgue  square  integrable.  At  first,  we  shall  restrict  our- 
selves to  functions  defined  over  a  finite  inteival,  which  for  simplicity  we  shall 
always  assume  to  be  the  interval  (0,1).  Our  space  p   is  then  the  space  of  all 
functions  u(x)  such  that 

J     u(x)  dx  <  00. 
o 

If  u  and  V  belong  to  /o  t   then  we  define 


*  "Bl  Friedman,  'Techniques  of  Applied  Mathematics  -  Theory  of  Distributions,' 
New  York  University,  Mathematics  Research  Group,  Research  Report  No.  EM-i|-7, 
October,  1952. 
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(3.2)  (u,^)  o /u(x)v(x)dx. 

o 


For  so.e  applications,  particularly  to  the  theory  of  Bessel  Unctions,   it  is 
comrenient  to  modify  the  definition  of  +v,«  JF  ^T>ions,   it  is 

i„^._.     .  aeiinition  of  the  space  5/  and  the  scalar  product  by 

introducing  a  non-negative  "weight*     function  ^(^)     J  ■ 
A-ctions  u(x)   such  that  ^^'  ^  ''  "°"  '''  ^^^^  ''  ^^ 


1 
y  u(x)  w(x)dx  <  00, 


and  the  scalar  product  becomes 

1 

(u,v)  =/u(x)v(x)w(x)dx. 
o 

^^  we  shall  be  Interested  in  linear  dlrre,.ntial  operators  s.ch  as  ^  or 

^  ^  *  S  -"ng  on  /.     These  operators  cannot  be  applied  to  e,ery  function 

if  ^    TT'  ""''  '""  '""'"°"  "°^  ""^  ^  dlri-er^ntlable  and,  second    eren 
Xf  the  f^ct^on  is  differentiable,  the  derivative  „a.  not  belong  ij 
let  L  be  the  following  linear  differential  operator: 

^°'=''S^^l''='£4*-.Vi(x)4.p„(x). 

»e  Shall  discuss  the  prTcise  .t^'oV;:  dTf  r:::::;.!:  ^te?^^™ 

wet,.  ^  ^^^^^^^_^  ^_^^^^^_^^  ^^  delta    :    10"    \et 

llToTr  f  '  =""  "^'  '^  -^-^  -  --"^=  -^-n  be  ;alL  It 

^.  f ;  ^"^iZtirtr: '" "-"- "  ^  ^^-^  --^^-^^  -^- 

.i„  .    .  ^PP-L^°^ti°ns,  the  domain  is  too  large  a  manifold,   since  it 

always  contains  an  infinite  n^ber  of  solutions  of  the  equation 


Lu  =  f(x). 


Put 
,„)        V">=gV^'(o,.gb^^„(3)(,, 

wHa«  u(J)(.,  represents  the  ,-th  derivative  of  „  if  ,  ,  ^  ^,  ^^^ 
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u^  '^(x)  «  u(x)  if  .1  =  0.  The  equations 

(3.3)  Bj^(u)  =0,   k  =  1,  2,  •••,  n 

will  be  n  linear  homogeneous  initial  or  boundary  conditions*  The  conditions 
are  linear  since  u  and  its  derivatives  appear  onlv  in  the  first  degree j  they 
are  homogeneous  because  if  u  and  its  derivatives  are  all  multiplied  by  a  con- 
stant c,  the  foiTO  B,  (u)  is  multiplied  by  the  same  constant  c.  The  conditions 
will  be  initial  conditions  if  all  b,  .  =»  Oj  they  will  be  end  conditions  if  all 
a.  .  =  Oj  and  otherwise,  we  will  call  them  boundsry  conditions.  For  our  present 
purposes  this  distinction  is  imnaterial.  We  must  assume,  however,  that  the 
conditions  (3.3)  are  independent,  that  is,  there  do  not  exist  constants  c^ ,  c^, 

•••,  c  such  that 
'  n 

CtBt(u)  +  c_B,(u)  +  ...  c  B  (u)  »  0 
11      t  ii  n  n 

for  all  functions  u(x). 

The  conditions  (3.3)  determine  a  linear  manifold,'^,  in  ^,  It  is  to 
be  noted  thatA^  is  not  a  closed  subspace,  that  is,  there  may  exist  in  />(_  a 
sequence  of  functions  u  (x)  which  converge  to  a  limit  u(x)  in  cj  ,  but  u(x)  is 
not  in  ??( ,  For  example,  let  TK.  be  defined  by  the  condition  u(0)  =  0.  Then 
consider  the  sequence  of  functions 

u  (x)  =0,  X  <  1/n 
=1,  X  >  1/n, 

It  is  clear  that  the  sequence  u  (x)  converges  in  the  mean  to  the  limit  function 
u(x)  =  1,  but  the  limit  function  does  not  itself  belong  to  '^. 

Logically,  we  should  consider  L  is  the  differential  operator  restricted  to 
that  part  of  its  domain  which  also  belongs  to  /?(_ ,  so  that  the  same  differential 
operator  with  different  sets  of  initial  or  boundary  conditions  represents  differ- 
ent linear  operators.  However,  for  convenience  in  writing,  we  shall  use  L  to 
represent  the  formal  differential  operator  only  and  use  L  and  the  appropriate 
manifold  to  represent  the  complete  linear  operator. 

It  will  occasionally  be  useful  to  consider  a  system  of  n  first-order 
differential  equations  instead  of  one  n-th  order  differential  equation.  In 
such  a  case  the  appropriate  space  is  the  space  of  vectors  U  with  components 
u-(x),  Up(x),«»»,u  (x),  where  each  component  is  a  real  square  integrable  function. 
The  appropriate  scalar  product  is 
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(V 


,U)   ^f   [u^v^*  U2V2+   •••   *^v^]dx. 


Note  that  this  is  a  ccMbination  of  the  scalar  product  in  an  n-dimensional  vector 
space  with  the  scalar  product  in  a  space  of  functions  of  one  variable. 

The  differential  operators  in  this  space  can  be  any  expression  of  the  form 

dx  dx 

where  now  A-, •••,!,  are  matrices  whose  elements  are  piecewise  continuous  functions 
of  X.  We  shall  usually  restrict  ourselves  to  the  consideration  of  the  simplest 
such  operator, namely 

consequently,  a  typical  differentiel  ec-uation  would  be 

(3.U  g  .  AU. 

Here  A  is  a  matrix  with  elements  a. .(x)  (i,j  =  l,2,»'»,n). 

When  equation  (3.U)  is  written  in  full,  it  looks  like  this: 

du^ 

dF=  ^1^1*  ^2^2  *  ••*  -^  ^n% 
•      •      •     •     • 

du 

■dF=  V^*  ^n2"2  ^  •••  *  ^nnV 

The  fact  that  we  consider  only  first-order  differential  operators  is  not  as 
restrictive  as  it  may  appear,   because  we  shall  now  show  that  every  n-th  order 
differential  equation 

u^"^  +  p^(x)u^"-^^   ...  +  pjx)u  =  0 

can  be  written  in  the  form  (3.U).  Put 

u^(x)  =  u(x),  UgCx)  =  u>(x),'**,  u^(x)  =  u^'^'^Vx) 

and  put 


no 


1 

0 


(3.5) 


A  - 


•Pn  -Pn-1 


0         1 
-P2     -Pi 


then  (3.1i)  is  equivalent  to  the  given  n-th  order  differential  equation. 

Similarly,  it  can  be  shown  that  a  system  of  n-th  order  differential 
equations  may  be  reduced  to  a  system  of  first-order  differential  equations. 

Of  course,  the  differential  operator  in  (3.U)  is  not  completely  defined 
until  we  specify  the  boundary  conditions  on  U.     Since  (3«l.i)  is  a  first-order 
system,  we  should  expect  one  vector  condition  on  the  values  of  U  at  x  ■  0 
and  X  =  1.     Since  one  vector  condition  implies  n  scalar  conditions  (one  for 
each  component),  we  may  write  the  conditions  on  U  as  follows: 


B^(U)  =  BgCU)  =   ...  B^(U)  -  0. 

Here,  B.(u)j,  (i  ■  l,2,«..,n)  is  a  linear  scalar  expression  containing  the 
values  of  u-,,  Up,  •'.,  u  at  x  -  0  and  x  ■  1. 

Adjoint  of  a  Differential  Operator;  Hermitian  Operators 

It  is  important  to  consider  simultaneously  a  linear  operator  ^  ^od  its 
adjoint  L  .  The  adjoint  is  defined  by  considering  the  equation 


(3.6) 

and  putting 


(v,Lu)  =  (w,u) 


w  ■  L  V, 


For  differential  operators,  an  equation  similar  to  (3.6)  can  be  obtained 
by  integration  by  parts.  For  example,  let  L  be  -s—  on  the  manifold'?^  de- 
fined by  the  condition  u(0)  »  2u(l).  Then 


Ill 


1 


du 


(v,Lu)  =y  V  ^  dx  =  vu  ^  -/u 


1    1 


(3.7) 


dx 

o 

1  *r 


dx 


u(l)  [v(l)-2v(0)]  -y^u  g  dx  -  (w,u) 


We  see  that  L  consists  of  two  parts:  a  differential  operator,  -  j~  5  and 
some  boundary  terms.  The  differential  operator  -  -r-  is  called  the  formal 
adjoint  to  the  differential  operator  -t—  •  The  adjoint  to  L  on  the  manifold 
9^  will  be  the  formal  adjoint  -  -r-  on  the  manifold  defined  by  the  condition 
v(l)  =  2v(0),  Now 


where 


(v,Lu)  =  (L  v,u) 


T»     dv 

L  V  »  -  •3— 
dx 


and  where  v  satisfies  the  condition  v(l)  =  2v(0), 

It  is  to  be  noted  that  in  this  example  L  acts  on  the  manifold  of  square 
integrable  functions  u(x)  such  that  u(0)  «  2u(l),  while  L  acts  on  the  mani- 
fold of  square  integrable  functions  v(x)  such  that  v(0)  ■  v(l)/2.  This  is 
an  illustration  of  the  general  situation  vdiere  the  manifold  on  which  L  acts 
ma^y  be  different  from  the  manifold  on  which  L  acts.  We  call  one  such  mani- 
fold the  dual  of  the  other. 

If  L  =  L  ,  the  differential  operator  is  said  to  be  formally  self-adjoint. 
If,  in  addition,  the  boundary  conditions  for  L  and  L  are  equivalent,  that 
is,  if  they  define  the  same  linear  manifold,  then  the  operator  is  said  to  be 
self-adjoint.  As  an  illustration,  suppose 

L  -  e^  1,  +  e^  f 

on  the  manifold  defined  by  the  conditions 

n'(0)  =  0,   u(l)  =  0. 

(v,Lu)  -/v(e\«»+  e\')  -  A vCe'Si' ) 'dx 
0  ^ 

1    1 
(3.8)  =  L^e^^'Jo  -y(v'e^)u'dx 

=  [ve^u']  -  [v'e^i]  +  /'u(v'e^)»dx 

«-  u'  (l)v(l)e^  +  V'  (O)u(O)  +  /■  uCe'^v"  +  e\'  )dx 


Then 
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This  shows  that 


In  order  that 


2 


(v,Lu)  -  (L  v,u) 


for  all  u,  the  boundary  conditions  that  t  satisfies  must  be 

v'(0)  »  0,  v(l)  =  0. 

Since  L  -  L  ,  the  differential  operator  is  foimaHy  self- ad  joint,   and  since 
the  boundary  conditions  on  u(x)  and  v(x)  are  the  same,  the  differential  oper- 
ator is  self-adjoint. 

The  adjoint  for  the  general  differential  operator  of  n-th  order 

.n  ,nvl 

L  =  Po^^^ -S  *  Pi^^^  TlJiT  *  ...  -^Pjx) 
dx  dx 

which  acts  on  the  manifold  defined  by  the  conditions 

Bj^(u)  -  0,  k  -  1,2,. -.jn, 

can  be  found  in  the  same  way.  Integration  by  parts  shows  that 

/"vLudx  «  /  ^JP  u   +  •••  p  u  J  dx 

*^a        a 

(3.9)  b 


yu[p^v-(p^_^v)'+  ...(-)^(p^v)^^^]dx  +  J(v,u)|^  . 


Here  J(v,u)  represents  the  terms  that  are  produced  during  the  integration 

by  parts.  For  example,  in  (3.7) 

.1 
J(v,u)l  »  u(l)v(l)-u(0)v(0) 

while  in  (3.8)       ^ 

J(v-,u)  ^  -  [ve^'  +  ve\]^  -  u(ve^)'l   . 

We  shall  call  J(u,t)  the  conjunct"^  of  the  functions  u  and  v.  It  is  clear  that 
J(u,v)  is  a  linear  homogeneous  function  of  u  and  v  and  their  deidvatives. 


+  Ince,  H.L,,  in  his  book,  Ordinary  Differential  Equations  (Dover,  19ltU), 
uses  the  term  "  bilinear  concomitant"  . 
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Equation  (3.9)  may  be  written  as  follows: 

1  1    ^  ,1 

/vLudx  *  /ul  vcbc  +  J(v,u)|      , 


where  L  v  is  the  differential  operator  defined  by  the  following  formula: 


(3.10)  l\  -  p^v  -   (Pj^_-^v)'  +   ...(-)"(p^v) 


(n) 


Now,  all  that  is  needed  to  complete  the  definition  of  the  adjoint  of  L  is 
to  determine  the  manifold  on  which  L  acts.  This  is  done  by  requiring 
JCuj'^)!  to  vanish  when  u  is  any  function  in  the  manifold  on  which  L  acts 
and  V  is  any  function  in  the  manifold  on  which  L  acts.  We  know  that  the 
manifold  on  which  L  acts  is  determined  by  the  boundary''  conditions  on  u, 
B,  (u)  «  0.  These  conditions  are  n  linear  equations  connecting  the  2n 
numbers  which  are  the  values  of  u(x)  and  its  first  (n-l)  derivatives  at 
X  »  0  and  x  =  1.  From  these  linear  equations, n  of  the  numbers,  say  u(0), 
u«(0),  •••  u^^^(O),  u(l),  u'(l),  •••  u^""^"^'(l),  can  be  expressed  in 
terms  of  the  n  mmbers  remaining,  namely  u     (O),  •••,  u  ~  (O), 
u^  ~  ~  (1),  •••  u  "  (1).  Substitute  these  expressions  in  J(v,u)|jj 
and  we  have 

u(^*l^(0)B*(v)  .   u(^*^^0)b:(v).  ....  u^"-^^(1)b;(v), 
i.  d  n 

where  B-(v),  B^Cv)**.**  B  (v)  are  n  linear  expressions  involving  v(x)  and 
its  first  (n-l)  derivatives  at  x  «  0  and  x  ■  1,  Consequently  the  mani- 
fold  on  which  L  acts  will  be  defined  by  the  boundary  conditions: 


Bj(v)  -  B2(v) B^(v)  -  0. 


The  examples  worked  out  at  the  beginning  of  this  section  are  illustrations 
of  this  method. 

In  the  case  of  a  system  of  first-order  differential  equations, the  ad- 
joint is  obtained  in  a  manner  similar  to  that  above.     Suppose  that  U  is  a 
column  vector  whose  n  components  are  functions,   and  supnose  that  A  is  an 
n  by  n  matrix  whose  elements  are  functions.     Consider  the  following  system 
of  first-order  differential  equations  : 

O.n)         §  .  Au 

with  the  boundary  conditions 

B,(U)  =  B,(U)  -  ...  B„(U)  -  0. 


nil 

Here,  B,  (U)  (k  =  l,2,»'*,n)  is  a  linear  scalar  expression  containing  the 
values  of  the  components  of  U  at  x  =  0  and  x  =  1» 

The  system  adjoint  to  (3,11)  is  obtained  by  taking  the  adjoint  of  tf.e 
operators  that  appear  in  (3.11) J  consequently,  if  we  suppose  V  is  a  coluinn 
vector  whose  n  components  are  functions,  then  the  adjoint  system  is 

(3.12)  -  S  -  A»V 

where  A     is  the  matrix  adjoint  to  A, 

The  adjoint  system  can  also  be  obtained  by  an  integration  by  parts, 
but  to  do  this  we  must  use  the  appropriate  scalar  product.     Suppose  that 
V   is  the  row  vector  whose  components  are  those  of  V,  and  suppose  that  V'»U 
represents  the  ordinary  scalar  product  of  vectors,  that  is,  the  sum  of  the 
products  of  corresponding  components}  then, as  we  have  defined  before, 


(V,U)   =  /  V    .  Udx 


Now,  put 


LU  = 

dU 
dx 

-  AU 

and  ' 

we  have 

1 
dx  -  y^V' 

(V,LU) 

dU 
dx 

•AUdx 

This 

shows  that 

the  adjoint  system  is 

d7' 
"  dx 

V«A  -  0. 

1^    >dV  } 

-  V'.U     '/    ~—  'Udx-  /  V'A.Udx, 
o/     dx  V 

'   ^0  ''o 


When  this  is  re-written  in  terms  of  column  vectors  instead  of  row  vectors  it 
becomes 


dx 


-  ^  -  A  V   =  0, 


which  we  have  already  obtained  previously  in  Equation  (3.12), 

The  boundary  conditions  on  V  are  found  from  the  condition  that  the  term 
which  appeared  in  the  integration  by  parts,  namely  V'»ul   ,  should  vanish 
identically  for  all  vectors  U  that  satisfy  the  boundary  conditions.  For 
example,  suppose  that  U  has  the  components  u^(x),  u«(x),  and  A  is  such  that 
equation  (3.11)  becomes 

suppose  also  that  the  boundary  conditions  are  as  follows: 
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(3.13)       u^(0)  -  Ugd),  UgCo)  -  0. 

Now,  let  V  have  the  components  v^(x),  Vp(x)j  then  the  equation  satisfied 
by  V  is 

dV 

To  obtain  the  boundsry  conditions  satisfied  by  V,  we  notice  that 

1  1 

V'.UI^  -  ^lV^2^2'o"  ^i(l)Ui(l)+V2(l)u2(l)-v^(0)u^(0)-V2(0)u2(0) 

(3.1U)  -  v^(l)u^(l)  +  [v2(l)-v^(0)]u-^(0), 

where  the  last  result  is  obtained  by  the  use  of  (3.13).  It  is  clear  that 
(3.1U)  will  vanish  identically  for  all  vectors  U  if  and  only  if 

v^(l)  -  0,  V2(l)  -  v^(0). 

These  are  the  required  boundary  conditions  on  V. 

The  discussion  of  thds  section  may  be  summarized  in  the  following 
Method.  To  find  the  adjoint  of  a  differential  operator  L  in  a  space  0 
consider  the  scalar  product  (v^Lu).  With  the  help  of  integration  by  parts, 
consider  it  as  the  scalar  product  of  u  with  some  vector  w,  which  depends  on 
y.  The  transformation  from  v  to  w  defines  the  adjoint  operator  L  .  The 
boundary  conditions  on  v  are  determined  by  the  requirement  that  the  terms 
resulting  from  the  integration  by  parts  vanish. 

It  is  to  be  noted  that  the  form  of  the  boundary  conditions  on  v  is  not 
unique;  for  if  P,(v)  =  pl^v)  ■  0  are  any  two  boundary  conditions  on  v  then 
the  conditions  a^p^(v)+a2P2('')  "  '^■3-M)-^^'^o^'^)   "  0  are  completely  equi- 
valent to  the  original  ones,  as  long  as  '^iPo'^pPl  ^  ^*     Nevertheless, 
despite  the  fact  that  the  boundary  conditions  are  not  unique,  it  is  easy 
to  show  that  the  manifold  defined  by  any  set  of  boundary  conditions  is 
always  the  same* 
Problems ; 

3.1  Find  the  adjoint  differential  operator  L  and  the  manifold  on  which 
it  acts  if 

(a)  Lu  -  if'+a(x)u'+b(x)u, where  u(0)  =  u'(l)  and  u(l)  =  u«(0), 

(b)  Lu  =  -(p(x)u')'+q(x)u  where  u(0)  =  u(l)  and  u'(0)  =  u'(l). 

Assume  that  the  scalar  product  is 
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1 
(u,v)  =  /^uv  dx, 
^o 

3.2  Suppose  that  Lu  =  u",  where  a,u(0)+b^u' (0)+c^u(l)+d^u' (l)  »  0  and 
a«u(0)+b2u'(0)+Cpu(l)+d2u'(l)  «  0.  Find  L  and  the  manifold  on  which  it 
acts.  For  what  values  of  the  constants  a.  ,b^, • ••^c^jd^  is  the  operator 
L  self-adjoint? 

3.3  Show  that  the  differential  operator 

Lu  -  -  -^^  (p(x)u')'  +  q(x)u 
is  formally  self-adjoint  if  the  scalar  product  is 

(u,v)  »  /  u(x)v(x)w(x)dx, 
3,U  Show  that  the  differential  operator 

Lu  =  a(x)u"  +  b(x)u'  +  c(x)u 

can  be  identified  with  the  operator  considered  in  Problem  3«3  by  putting 
q  =  c,  p  -  exp  /  b/a  dx,  w  = -e  exp  | -/ b/a  dxj  ,  It  follows  that  the  general 
second-order  differential  operator  is  formally  self-adjoint  if  the  appro- 
priate scalar  product  is  used, 

3.5  Find  the  system  adjoint  to  the  following; 

Pi   P2   0 

^1       ^2       ° 
where  u^(0)+U2(0)+u-(0)   ■=  0,  u-|^(l)   =  UgCD  and  u^(l)  »  u^(l)+U2(l). 

3.6  Consider  the  complex-type  scalar  product 

1  

(v,u)  =  f  v(x)  u(x)w(x)dx, 

where  w(x)   is  a  non-negative  real  function.     Is  Lu  -  -  -  (pu')'+qu  formally 
self-adjoint  in  a  space  with  the  above  scalar  product?     Distinguish  between 
the  case  where  p  and  q  are  both  real  functions  and  the  case  where  p  and  q 
are  complex  functions. 
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l«lf-Ad,1o)nt  Sec<>nd-0r.1.r  ln^fer.ntlal  np.,.*... 

seif-adjoint  if  the  appropriate  scalar  product  is  used-  ty,«.  ^ 
we  ahnii  w■..^^■^  iu  v/u.w,i,  xb  usea,  tnerefore 

we  snail  write  the  operator  as  follows: 

Lu  -  -  -  (pu ' )  •  +  qu, 

case.  ::::: '" :"  "-•  "^ '-  ">•  -"---  -  -  - '-« -  =er.. 

ope^tor  L  1.   ealu.  pcltWe-deflaUe   If   („.^)  >  o  .„  all  ,al.e.  of  „  except 
«  -   0.     Nov,  .h.„  Z  1,  defined  V   (,.15),  „,  have 

(u.ln)  -y^C-l  (pu')'  +  q^   ,4, 
0 

=y    (pu'^+  qwu^)   dr  -  puu'    I      ; 

consequently.  If  p  >  0  ,  q  >  o  w  >  0  An^  <*•  *u  v 

o„^v  ,,  ^    ^         .  ^  -^  u.  w  >  0,  and  If  the  boundary  conditions  on  u  are 

..  JLuI'   ',r  n''  '°™"''  "^'-^J°'"*-      *^"  «^«  -Edition,  .m  t 

are  l^o.^a.'o.  ap:^:!:  '"''  "'  '"''^^  "^^  ""'""  '""  -'"  ""- 

1)      The  boundary  conditions  are  unmixed   If  ««.},  >.       ^ 
the  «1..,  Of  u  and  Ue  de.l,atl,e.  ,^^^7."  „     ^     "  iT'"""  '""'"' 

%^(o)  +  e„u'(o)  =  0. 
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It  is  easy  to  show  that  if  u  satisfies  an  unmixed  boundary  condition  at  x  =  0 
and  an  unmixed  boundaiy  condition  at  x  =  1,  then  L  is  self-adjoint. 
2)  The  boundary  conditions  are  periodic  if  they  have  the  form: 

u(0)  =  u(l),  u'(0)   =  u'(l). 
In  this  case  again  it  is  easy  to  show  that  L  is   self -adjoint . 
We  now  prove 
Theorem  I .     If  u  is  any  solution  of  the  equation  Lu  =  0  and  if  v  is  any  solu- 
tion  of  the  equation  L  v  =  0,   then  the   conjunct  of  u  and  v  is  a  constant  whose 
value  depends  on  u  and  v. 

/Here  L  and  L     are  formally  adjoint  differential  operators  in  the  full  space 
The  proof  is  an  immediate  consequence  of  the  definition  of  the  conjunct. 

From  equation  (3.9)  we  have 

b  b      _^  b 

/'vLu  dx  -  /*   uL     V  dx  =  J(u,v)    | 

3i  3.  at 

From  the  hypothesis  on  u  and  v  it  follows  'ohat 

b 
J(v,u)    I     =0; 
a 

consequently   the  value  of  j(v,u)   at  x  =  b  is  the   same    as  its  value  for  x  =  a. 

Since   a  and  b  are  arbitrary  points,   this  proves  that  the  value  of  3(v,u)  is 

constant . 

Corollary.     If  L  is  a  formally  self-adjoint  operator  and  if  u,    and  Up  are  any 

solutions  of  Lu  =  0,   then  the  conjunct  of  U-,    and  Up  is  a  constant  whose  value 

depends  on  the  functions  u-.    and  Ug    . 

Note  that  the  theorem  and  corollary  are  also  valid  for  an  n       order  system. 

As  an  illustration  of  this  theorem  consider  the  following  differential 
equation: 

(3.17)  u"  +  u  =  0. 

Since    (3.17)  is  fonnally  self-adjoint,  v  will  also  be  a  solution  of  it.     The 
conjunct  of  v  and  u  is  now 

J(v,u)  =  vu'   -  uv'    . 
Theorem  I  states  that  this  expression  is  constant  for  any  two  solutions  of  (3.16). 
For  example,   if  u  =  sinx  and  v  =  cosx,then 

J(v,u)  =  cosx  cosx  -  sinx(-sinx)  =  constant. 

This  is  obviously  correct  because 

2  .    2        - 

cos  X  +  sm  X  =  1 

for  all  values  of  x. 
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The  V.^ronsklen  of  the  Solutions  of  a  Differential  Equation 

Equ&tion  (3.16)  shows  that  for  the  self-edjolnt  differential  operator 
defined  in  (3.3.5)  the  conjvmct  Is 


(3.18) 


J(v,u)  =  p(x)(vu'-nv'). 


The  expression  in  parentheses,  vu'-uv',  is  called  the  Wronskian  of  v  and  u  and 
plays  8  very  important  role  in  the  theory  of  differential  equations.  We  shall 
denote  the  Wronskian  of  v  and  u  "by  the  sym'hol  [j,-q\  .  The  following  properties 
of  the  Wronskain  are  evident: 

(a)      [v,u]   =  -   [n,v]  ; 

(t)  for  any  constants  ou.   and  ag.  we  have 

(c)  for  any  three  functions  u(x),  v(x)  and  w(z)  we  have 

u[y,w]   +  v[w,u]   +  w[u,v]    =    a. 

The  last  property  follows   immediately  from  the  expansion  of  the  following  deter- 
minant : 


u 
u 
u' 


V 
V 


w 
w 


=  0 


The  following'  theorem  shows  why  the  Wronskien  is  important: 
Theorem  II.  If  the  Wronskian  of  two  ftmctions.  u(x).  v(i)  is  identically  zero. 
then  the  functions  are  linearly  dependent   in  any  intervel  where  they  do  not 
vonish;conver8ely.  if  the  functions  are  linearly  dependent,  their  Wronskian  is 
identically  zero.  V^e  prove  the  converse  first.   Suppose  there  exist  constants 
or.,  P,  not  "both  zero,  such  that 


(M9) 


rTu(x)  +  Pv(x)  =  0 


for  all  X  in  an  Interval  I.  Then,  dlfferentiRting  (3.19)  with  respect  to  x, 
we  have 


(-^.20) 


rni'(x)  +  ew'(x)  =  0 


for  all   X  in   I.     Iquntions    (>.19)  and    (3.20)   may  be  considered  as  homogeneous 
equations  for  the  quantities  a  and  P.      Since  the  equations  have  a  non-triviel 
solution,    the  deterrninsnt   of  the  coefficients  must  be   zero,   that  is, 
u(x)  v(x) 


u'(x) 


(x) 


=   [u.v]      =  0 
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This  proves  that  the  Wronsicien  is  zero  for  all  x  In  I. 

To  prove  the  first  part  of  the  theorem,  suppose  that 

Ei.-^  =  ■a(x)v'(T)-v(x)u'(x)  =  0 
for  all  r   in  I.  Then,  if  v(x)  and  n(x)  are  not  zero  in  I,  we  have 


v_ 

v( 


which  can  "be  Integrated  to  give 

log  v(x)  =  log  u(x)  +  log  c 
or 

v(t)  -   cu(x)  =  0 

for  all  X  in  I.  This  shows  that  v(x)  and  uCx)  are  linearly  dependent. 

From  Theorem  I  and  (3.18)  it  follows  that  if  u^  and  Uj,  are  any  solutions 
of  the  equation 

L  ti  =  0 
where  L  is  defined  "by  (3.15),  then 

(3.21)  ^ti^.ti^  =  u^u^  -  UgtLj  =  const.  /p(i) 

for  0  <  X  ^  1,  The  value  of  the  constant  will  depend  upon  the  particular  solu- 
tions tL.(i),  U2(x)  that  are  considered. 

Notice  that  in  (3.2l)  we  have  implicitly  assuaed  that  p(x)  4  0   for 
0  <  X  <  1.  A  point  X  such  that  p(x  )  =  0  is  called  a  singular  point  of  the 

0  0  ~ 

differential  equation.  Let  us  re-write  (3.2l)  in  the  following  form: 

(3.22)  p(x)(u^u^  -  '^2^V  ~   ^°^"^*' 

The  vnlue  of  the  constant  can  ^e  determined  for  any  value  of  x  ^f  x  .  If  u^  and 
u  are  not  linearly  dependent,  the  value  of  this  constant  will  not  he  zero  unless 
p(x)  =  0.  I'That  happens  to  (3.22)  as  x  approaches  x  7   Since  p(x)  approaches  zero, 
while  the  product  remains  a  non-zero  constant,  we  conclude  that  at  least  one  of 
the  functions  •ul,  (x)  and  u„(x)  must  hecome  infinite  as  x  approaches  x  .  later 
we  shall  disff>iss  in  great  detail  the  nature  of  this  singular  behavior  of  the 
solutions  of  the  differential  eqiiation. 

The  result  (3.2l)  has  heen  derived  only  for  solutions  of  a  second-order 
differential  equation.  There  is  an  analogous  result  for  solutions  of  an  n-th 
order  differential  equation  or  for  solutions  of  a  system  of  differential  equations. 
Before  stating  the  result,  we  must  first  extend  the  definition  of  the  Wronskain 
80  that  it  will  apply  to  a  set  of  n  fiinctions.  The  Wronskian 
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of  ii,(x),  UgCx),  ...,  u^(x)  is 


(3.23)  W(u^,U2,...u^)  = 


^     ^2 


u 


u. 


n 

I 

1 
n 


\ 


(n-1)  „(n-l) 


u; 


u 


(n-1) 


n 


'i'he  primes  and  the  superscripts  denote  differentiation  with  respect  to  x. 

We  now  state 
Theorem  III.  If  u, ,Up  ...,u  are  arbitrary  solutions  of  the  differential 

equation 


(n)      (n-1) 
u^  '  +  p_u^ 


+  ...  + 


p  u  =  0. 
^n     * 


then 


^p^dx 
W(u-^,U2,...,u^)i:c  e      , 


where  c  is  a  constant  depending  on  t he  particular  solutions  u, ,Up , . . . , u  . 

This  theorem  will  be  an  immediate  consequence  of 
Theorem  IV.  If  U-, ,  U^,  . . .  U  are  arbitrary  vector  solutions  of  the  system 
of  n  differential  equations 

dU 


(3.2U) 


dx 


AU 


and  if  we  define  the  Wronskian  W  as  the  vaLue  of  the  determinant  whose  columns 


are  the  vectors  U, ,  U^    ...»  U   , 
1'     2         '     n 


then 


+/trace  Adx 


W  =  ce 
where  c  is  a  constant  whose  value  depends  on  the  particular   solutions 

"1*^2    •••  ^n   • 

In  order  to  avoid  the  difficulties  of  a  cumbersome  notation  we  shall 
prove  this  theorem  for  the  case  n  =  2.     The  general  case  can  be   treated  in  a 
similar  way. 

SupTDO  se   that 


and  that 


A  = 


a      b 


c       d 
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then 

We  hare 
(■-^.25) 


W  = 


^1       ^2 


w         w 

1  2 


dW 
dx 


w. 


tJsing   (3.24)  we  ohtain 


w 


and  a  corresponding  result  for  v.,  w^.     When  these  results  are  used  in   (3.25) 
we  get 


*1 

dx 


11  2        2 


cv^+dw^        cv -+dw 


=   (a+d) 


w. 


=  W  trf;ce  A, 


and  then  an  integration  gives  the  conclusion  of  Theorem  17: 


w=  C 


/. 


traceAdr 


Theorem  III  follows  from  Theorem  IV  if  A  la  defined  as  in  (3.5). 
There  is  an  interesting  consequence  of  TheoremsIII  and  IV  which  we  state  as  a 
Corollary,   Suppo-^e  that  in  Theorem  IV  trace  A  is  integrable  over  the  interval 
(s-.h);  then  if  the  Wronalcian  vanishes  at  one  point  of  the  interval  it  vanisbes 
identically  in  the  whole  interval  and  J;he  vectors  IL,,*-',U  are  linearly  dependent. 
The  proof  is  easy,  "Prom  Theorem  IV  we  know  that 


/ 


trAdx 


If  W  vanishes  anywhere,  C-  must  "be  zero  since  the  exponentials  never  vanish; 
consequently  W  vanishes  identically.  Since  W  vanishes  for  all  x  in  (a,h),  there 


exist  functions  a-(x). 


,  a  (x)  such  that 
n 
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a,U,  +   •••  +  a  U     =  0   , 
11  n  n 

for  all  X  in   (e.^).     At  least  one  of  the  a   (x),   say  a   (x),    is  not   identically  lero; 
therefore  we  ipay  write 

^1  =  V2-^    "'*\\' 

where  K  =  -«.'»,.  (k  =  2,*'*,n).  Substituting  this  expression  for  U  into  (3.2J^) 
we  fet 

Either  all  the  derivatives  are  zero,  in  which  case  h  , • • 'b  are  constants  and  so  U^ 

^     °     dbp  ^ 

is  linearly  dependent  on  !!-,•• ',11  ;  or  at  least  one  of  them,  say  z. ,  is  not  zero, 

2     ^  dx 

in  which  c^se  we  may  write 

db,    dbj, 

where  c.  =  -  -r-^  I  -r^   •  (k  =  ^.•••,n}.  Now  we  t)roceed  as  above  and  we  show  that 
k     dx  '  dx 

either  C-,«'«,c  are  constants,  in  which  case  U  is  linearly  dependent  on  U^, . . . ,U  , 

or  U_  can  be  expressed  in  terms  of  U.,*'*,U  .   It  is  clesr  that  this  reduction  can 

be  continued  until  we  prove  the  corollary. 

Problems  t 

3.7     Prove   that  Bessel's  differential  operator  defined  by 

Lu=  -  ^  (xu')     +  "2 
X 

is  positive  definite  if  the  boundary  conditions  are  u(o)  =  0  and  u(l)+  hu(l)  =  0, 
where  h  is  any  positive  constant.   Hint:   Use  w(x)  =  x  as  a  weight  function. 

3.8.  Prove  that  the  operator  defined  in  (3.15)  is  self-adjoint  if  the  boundary 

conditions  are  unmixed  or  if  the  boundary  conditions  ere  periodic. 

(iv'l    2 

3.9.  Verify  Theorpn  I  if  Ln  =  u   '  +  k  u. 

3.10.  Verify  Theorem  I  if  tj^  is  the  space  of  two-component  vector  functions  U,and 


if 


^-  s  -  Co  :]  «• 


where  a,  b,  c  are  cons  ♦-.ants. 
3.11.   Prove  that 
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2         2 

3.12  Consider  the  fimctions  u(x)  «=  x  ,  v(x)  =  x  agn  x  over  the  interval  (-1,1), 

The  Wronskian  of  u  and  v  is  identically  zero,  yet  u  and  v  are  not  linearly  depen- 
dent, Explain  the  apparent  contradiction  to  Theorem  II.  Hint:  u  and  v  vanish 
when  X  =  0. 

3.13  Stippose  that  La  =  -Cxa  )  .  What  is  the  Wronskian  of  any  two  solutions  of 
Iai  =  0?  Investigate  the  "behavior  of  the  Wronskian  as  x  approaches  zero. 

3.14  Prove  Theorem  III  directly.  Hint:  Use  a  method  similar  to  that  used  in 
the  proof  of  Theorem  IV. 

Green's  Function 

The  inverse  of  an  operator  L  Is  defined  to  he  operator  L"  such  that 

L'-'-L  =  I, 

I  heing  the  identity  operator.  In  E  ,  I  is  represented  by  a  matrix  with  elements 
6   which  are  defined  as 

*ij  =  °'     ^     *     ^ 
=  1,      i  =  J. 

In  /3  ,  I  will  he  represented  not  hy  a  me.trixi  "but  "by  an  integral  operator  whose 
kernel  is  a  symbolic  function,  namely  6(xf-t),  which  equals  zero  when  x  j^  t  and 
which  is  -undefined  in  value  when  x  =  t.  Therefore,  in  the  case  of  linear  differential 
operators,  the  inverse  will  he  found  "by  solving  the  differential  equation 

(3.26)  Ig  =  6(x-t), 

subject  to  the  given  boundary  conditions.  Here,  g  is  a  symbolic  function  for  the 
distribution  in  x  and  t  which  satisfies  (3.26).  It  will  be  called  the  Green's 
function  of  the  differential  operator  L  in  the  domain  defined  "by  the  given  boundary 
conditions. 

Since  6(x-t)  is  zero  in  the  neighborhood  of  every  value  of  x  except  x  =t  , 

(3.26)  implies  that  g(x,t)  satisfies  the  homogeneoias  equation 

(3.27)  l€  =  0 

except  when  x  =t  .  In  the  neighborhood  of  x  =  t,  the  behavior  of  g(x,t)  will 
depend  on  the  explicit  form  of  L, 

Suppose  L  is  a  second-order  self-adjoint  operator: 

(3.15)  1  =  -  aJ  (P  i)  ♦  \u). 

and  assume  p(t)  /  0  in  (3.26).  Put 

(3.28)  g(x,t)  =  -(x-t)H(x-t)/p(t)  +  k(x,t). 
Now, 
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dx 


Q(x"t)H(x-t)  ]   =  H(x-t) 


and  , 

2-  H(x-t)  =  6(x-t)   ; 

consequently  from  (3.28)   and  (3.l5),  we  have 

(3.29)  Lk   =  q(x)(x-t)H(x-t)/p(t)    . 

Note  that  the  right  side  of  (3.29)  is  a  continuous  function  of  x.  It  can  be 
shown  that  k(x,t)  will  be  a  twice  differentiable  f\inction  for  all  values  of  x 
such  that  p(x)  ^0.  In  fact,  we  shall  prove  the  following  theorem  in  Appendix 
II: 

Theorem  II-I .  Let  p(x),  q(x),  and  f(x)  be  piecewise  continuous  functions  of  x 
in  the  closed  interval  (0,l),  and  let  p(x)  be  positive  in  that  interval;  then 
there  exists  a  continuous  function  u(x)  such  that  pu'  exists  and  is  continuous 
for  all  x;  furthermore 

(II. I)  (pu')  -  q(u)  =  f(x) 

for  all  vaJ.ues  of  x  for  vrhich  both  sides  are  continuous  functions  of  x;      and 

u(0)   =  u'(0)  =  0. 

Since  k(x,t)  is  tvrice  differentiable,  equation  (3.28)   enables  us  to 
state  the  following  theorem  about  the  behavior  of  the  Green's  function: 

Theorem  V :  For  all  vaJues  of  x,  except  x  =  t,    the  Green's  function  satisfies  the 
homogeneous  equation.     In  the  case  where  L  is  given  by  (3.15),  the  Green's  func- 
tion at  X  =  t  is  continuous  while  its  derivative  has  a  jump  of  magnitude  -l/p(t), 
where  p(t)   f  0. 

Once  the   Green's  function  is  known,  the  non -homogeneous  equation 

(3.30)  Lu  =  f(x) 

with  the  assigned  boundary  conditions  can  be  solved.  The  solution  is 

1 

(3.31)  u(x)  =  /  f(t)g(x,t)  dt. 

0 
The  proof  that  u  as  defined  in  (3.31)   satisfies  (3»30)  is  as   follows: 

1  1 

Lu  =y^   f(t)Lgdt  =»  /*f(t)5(x-t)dt  =  f(x). 

u  also  satisfies  the   assigned  bo-undary  c  onditions  because  g(x,t)   as  a  function 
of  X  satisfies  them. 
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Example  One  -  Green's  7-anctlon 

dx 
and  let  the  "boimdary  conditions  "be 

ti(0)  =  u'Co)  =  0, 

To  find  the  Green's  fiinctlon  g(x,t)  we  must  solve 

(3.32)  i-f«.6(x-t) 
dx 

with  the  conditions 

(3.33)  g(0)  =  g^io)   =  0. 

As  remarked  prevlotisly,  g  Is  a  solution  of  the  homogeneous  equation 

(3.3^)       H=0' 

except  at  x  =  t,  where  It  la  continuous  and  its  derivative  has  a  jump  of  mag- 
nitude -1.  For_x_<^t,  (3.33)  and  (3.3^)  imply  that 

while  for  x  >  t  we  assume  that  g  is  an  arbitrary  solution  of  (3.3U),  namely 

e  =  cue  +  $, 

Since  the  .lump  in  the  derifative  of  g  is  -1,  we  mast  have 

a  =  -1; 
since  g  is  continuous  at  x  =  t,  we  must  have 

P  =  t. 

Tlnally 

g  =  0,   X  <  t 
=  t-x,  X  >  t 
or 

(3.35)       g(x.t)  =  -(x-t)H(x-t) 
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for  all  values  of  x.   It  is  oIdvIous  that  g(x,t)  satisfies  (3.32)  and  (3.33), 
Now,  the  solution  of 

(3.36)  -  ^  =  f  U) 
with  the  conditions 

(3.37)  u(0)  =  u'(0)  =  0 

may  he  fotmd.  Using  (3.31)  «^nd  (3.35),  we  have 

1 

(3.38)  ti(x)  =  -  /  f(t)(x-t)H(x-t)dt 

X 

=  -/  f(t)(x-t)dt, 
0 

and  this  is  clearly  the  solution  of  (3.36)  satisfying  (3.37). 

Suppose  that,  instead  of  (3.17) »  the  conditions  are 

u(0)  =  P,  u'(0)  =  a. 

Then,  put 

u  =  cat  +  P  +  v; 

we  find  that 

dx 
where 

v(0)  =  v'(0)  =  0. 

Solving  this  "by  means  of  (3.38)  we  find  that 

X 

u(x)  =  ax  +   P  -  /  f (t)(x-t)dt. 
0 

Example  Two  -  Green's  Function 

^2 

Let  L  again  he  -  — r  ,  hut  now  suppose  the  conditions  are 
dx 

(3.39)  n(0)  =  0,    u(l)  =  0. 

To  find    the  C-reen's   function  we  inust   solve 

^=-   6(x-t) 
dx2 


-  128  - 

on  the  manifold  defined  "by   the  hotindary  conditions 

f(0,t)  =  ^(l,t)  =  0. 

We  flhflll  illustrate  a  very  nseftil  technique  for  finding  the  Green's  function, 
^or  I  <"  t,  the  solution  of  the  homogeneous  equation  that  satisfies  the  left 
boundary  condition  g(0,t)  =  0  will  he  proportional  to  x.  for  x  >  t,  the  solu- 
tion of  the  homogeneous  eque.tion  that  satisfies  the  right  boundary  condition,. 
g(l,t)  =  0,  will  he  proportional  to  1  -x. 
Write 

g(x,t)  =  X  ,    X  <  t 
=  1-x,    X  >  t. 

This  is  wrone:  because  g(x,t)  is  not  continuous  for  x  =  t.  The  value  of  g  as 
X  approaches  t  from  below  is  t,  while  the  value  of  g  as  x  approaches  t  from  above 
is  1  -  t.  Multiply  the  first  expression  for  g  by  the  value  of  the  second  eaipression 
at  X  =  t  and  multiply  the  second  expression  for  g  by  the  value  of  the  first 
expression  at  x  =  t.  Then  write 

(3.^0)        g(x,t)  =  x(l-t),   X  <  t 

=  (l-x)t,   X  >  t. 

This  function  is  continuous  at  x  =t.  Because  the  derivative  of  this  function  for 
X  <r  t  Ir  (l-t)  while  the  derivative  for  x  >  t  is  -t,  the  Jump  in  the  derivative 
at  X  =  t  is 

_  t  -  (l-t)  =  -1  . 

(''.^O)  is  therefore  the  Green's  function  for  the  operator  -  — x  in  the  manifold 

dx 

defined  by  the  conditions  (3.39).  We  oey  write  the  expression  for  g(x,t)  as 

follows : 

g(x,t)  =  x(l-t)H(t-x)  +  (l-x)tH(x-t) 

-   2   -  xt  -   2 

Note  that  g(x,t)  is  symmetric  in  x  and  t.  We  shall  see  later  that  this  is 
generally  true  for  the  Green's  function  of  a  self -adjoint  operator. 

By  differentiation  we  may  check  that  (3.iJ0)  is  actually  the  solution  of  the 
dii'ferential  equation.  We  have 
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^=  -  6(x-t). 

As  in  the  preceding  example  we  find  that  the   solution  of 

-  ^    u  =  f  (x) 

with  the  ■boundsry  conditions 

n(0)  =  ud)  =  0 


is 


u  =  Af(t)[x(l-$H(t-x)  +   (l-x)t  TT(x-tldt 
0 

=   {\-T)ft  f(t)dt  +  xy    (l-t)f(t)dt. 


0 

Example  Three  -  Green's  Function 
Consider  the  operator 

dx'^ 

in  the  manifold  defined  "by 

u(0)  =  u'(l)  =  0. 

The  Qreen's  function  is  found  "by  solving  the  equation 

,2 

(?.iH)  S^+  iig  =  -  6(x-t), 

dx 

where 

(^.^2)  p(0,t)  -  g^(l.t)  =  0. 

We  use  the  same  technique  as  in  Example  Two.   The  solutions  of  the  homo- 
geneous equation  corresponding  to  (3.i<l)  are  any  linenr  com"bination  of  sin  2x 
and  cos  2x.  We  pick  two  solutions  of  the  homogeneous  equation  -  one  to  satisfy 
the  "bonndary  condition  at  zero,  the  other  to  satisfy  the  boundary  condition  at  unity, 
We  may  then  write 

g(x,t)  =  sin  2x       X  <  t 
=  cos  2(l-x),   X  >  t. 

Since  g(x,t)  is  not  continuous  for  i  =  t,  we  multiply  the  first  expression  "by 
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the  value  of  the  second  expression  at  t  =  t  and  also  multiply  the  second  ex- 
pression "by  the  value  of  the  first  at  x  =  t.  We  get 

p(x,t)  =  sin  2x   cos  2(l-t),  x  <  t 
=  cos  2(l-x)  sin  2t,   x  >  t. 

Now,  g(x,t)  is  continuous,  hut  the  Jump  in  the  magnitude  of  the  derivative  at 
X  =  t  is 

2  sin2(l-t)  sin  2t  -  2cos  2t  cos  2(l-t) 

=  -2  cos  2(t+l-t)  =  -2  cos  2. 

Since  it  follcws  from  (3,^l)  that  the  jump  in  the  derivation  should  he  -1,  we  see 
that  this  last  formula  for  e(x,t)  is  still  wrong.  To  correct  it,  we  must  divide 
it  hy  2  cos  2,  We  have  then  the  final  and  correct  results: 

g(x,t)  =  sin  2x  cos  2(l-t)/2  cos  2,   x  <  t 
=  cos  2  (l-x)  sjn  2t/2  cos  2,  x  >  t 

or 

(^.UJ)         ^(^^)^  sin  2x  cos  2(l-t)   ^^^_^^ 

c   cos  Z 

^  cos.^(l:od^in_2t  jj(   J 
2  cos  2 

Note  fi^ain  that  g(x,t)  is  symmetric  irt  x  and  t. 

It  is  easy  to  show  hy  straightforward  differentiation  that  g(x,t)  as  given 
iy  (3.^3)  satisfies  (3>l).  We  have 

dg  _  2  cos  2x  COS  2(l-t)  „/._  \  _j_  2  sin  2(l-x)  sin  2t  „/    s 
dx       2  cos  2  ~^  2  cos  2        ~ 

sin  2x  cos  2(l-t)  -  / .   v  .  cos  2(l-x)  sin  2t  ./  .  v 
-  2  cos  2     Mt-x)  +      2r^^^  6(^-t). 

The  sum  of  the  last  two  terms  is  zero,  hecause 
f(x)6(x-t)  =  f(t)6(x-t) 


and 

sin  2x  cos  2(l-t)  j     _  cos  2(l-x)  sin  2t 

x=t 


sin  2x  cos   2(l-t)    j  _     cos   2(l-x)   sin  2t 

2  cos   2  x=t     ~  2  cos   2 


Differentiating  the  formula  for  7^  ,  we  find  that 

dx 
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^  =  -Zl  sin  Zx   cos  2(l-t)  ^^^^^  ^  ^  cos  2(l-z)  sin  2t  ^^^_^^ 

.    c.  £,   COS    <;:  c   COS    c 

dz 

cos  2x  COS  2(l-t)      .  / .      V    .    sin  2(l-x)   ain  2t   -/      .\ 

=     _  %  +  Sls^tl   j-g^j^  2(l-x)  Sin  2t  -  COS  2x  cos  2(l-t3 

=     _  iig  -  6(x-t), 

which  proves  that  g(x,t)  satisfies  (3.^l). 

Green's  F-gnction  in  the  Case  of  Unmixed  Boundary  Conditions 

The  technique  that  has  been  lased  in  the  two  previous  examples  can  be  applied 
to  any  second-order  differential  operator  if  the  "boundary  conditions  are  not 
mixed,  that  is,  if  each  boundary  condtion  involves  just  one  point  of  the  boundary 
but  not  both.  A  condition  such  as  u  (o)  =  u(o)  is  not  mixed  but  a  condition 
such  as  u(l)  =  u(0)  is  mixed. 

We  shall  obtain  the  Green's  function  for  the  general  self-adjoint  operator 

1  =  -  A  (p  A) . , 

when  the  domain  is  defined  by  the  general  unmixed  boundary  conditions 

(3.Wj)  B^(u)  =  B^Cu)  =  0. 

We 'shall  assume  p(x)  ^  0  in  the  interval  (0,l)  and  that  B  (u)  involves  values  of 
u  and  its  derivatives  at  x  =  0  only,  while  Bp^^^^  involves  values  of  u  and  its 
derivatives  at  x  =  1  only. 

The  Green's  function  is  the  solution  of 

(3.^5)  (pg^)^  -  qgU.t)  =  -5(x-t) 

with  the  boundary  conditions 

(3.^6)  B^(g)  =  B^Cg)  =  0. 

Consider  the  homogenecos  equation 

(3.^7)  ^P^x^x  -  <ri  =  0 

without  any  boundary  condtions.  S'j^jpose  that  v  (x)  and  v^Cx)  are  any  two  inde- 
pendent sol-ations  of  (3.^7).  Let  w  (x)  be  a  linear  combination  of  v^  (x)  and 
v_(x)  which  satisfies  the  condition 

B^(w^)  =  0. 


i 
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Similarly,  let  vAx)   "be  a  linear  coa^bination  of  v^(x)  and  v^Cx)  which 
satisfied  the  condition 

BgCw^)  =  0  . 

Then  we  start  the  construction  of  the  Crreen's  fiinctlon  "by     writing 

g(x,t)  =  w^(x),   T  <  t 

=   WgCx),    X  >  t    . 

This  expression  for  g(x,t)  will  satisfy  the  differential  equation  (j.h^)   for 
X  ^  t  and  also  the  "boundary  conditions  (3.^6)^  hut  it  will  not  satisfy  the 
continuity  or  the  ^ump  conditions  at  x  =  t.  To  make  g(x,t)  continuous,  multiply 
the  first  expression  "by  the  value  of  the  second  at  x  =  t  and  multiply  the  second 
expression  by  the  value  of  the  first  at  x  =  t  so  thst  we  h-^ve 

(3.U8)        g(x,t)  =  w^(x)w2(t),  x<t 

=  w  (x)w.  (t),   X  >  t  . 

This  expression  still  does  not  satisfy  the  jump  condition  in  the  derivative  at 
X  =  t.  The  Jimp  in  the  derivative  at  x  =  t  is  given  "by 

w^(t)wj^(t)  -  w^(t)w2(t)  =  [w^.wj 

X  =  t 

while  it  should  he  -l/p(t).  This  can  he  remedied  hy  dividing  (3.^)  "by  the 
conjunct  of  Wj,,  w  ,  The  final  and  correct  formnila  is  then 

(3.^9)        g(x,t)  =  w^(x)w2(t)/|w2.wjp(t),   x<t 

-  W2(x)w^(t)/^2»^llp(*)»   3c>t, 

or,  written  in  another  way, 

w, (x)w.(t)H(t-x)  +  w-(x)w, (t)H(x-t) 
(3.50)  g(x.t)  =  -^ ^ -2 i 


[wg.Wjl  p(t) 


The  conjunct  in  the  denominator  of  these  expressions  is  to  "be  evaluated 
at  X  =  t.   However,  as  we  ha-ve  shown  in  (3.22) 


p(x)[w2»*'t1  =  constant; 


consequently  instead  of  evaluating  the  conjunct  at  x  =  t ,  we  may  evaluate  it  at 
any  convenient  point.  "P'or  a  simple  illustration  of  the  usefulness  of  this  fact, 
consider  Example  Three  of  the  Green's  function.  There,  the  conjunct  of 
cos  2(1-x)  and  sin  2x  had  to  he  found.  Now,  since  p(x)  =  1,  we  have 
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[cos  2(l-x),   sin  2x\   =  constant 

=  2  cos  2(l-x)  COB  2i  +  2  sin  2(l-x)  sin  23: 
=  2  cos  2 

when  we  put  x  =  0.     Before,   this  reevilt  was  olstalned  by  using  the  addition 
theorem  for  the  cosine,  hut  now  it  is  an  immediate  consequence  of  the  proper- 
ties of  the  conjunct. 

Ponnulas    (3.^9)  and   (3. 50)  break  down  if   [wo'^'il   ~  ^*     "^^  ^^^*  *^*   ^® 
conjunct  is   zero  implies   that  one  function,  say  vAx)t   is  a  multiple  of  the 
other,  w   (z).     Since  B^(w^)  =  0  and  B^Cwg)  =  0,    it  follows  that  B^Cw^) 
=  Bp(wp)  =  0.     This  means  w^Cx)  is  non-trivial  solution  of  the  homogeneous 

equation  i 

(pu  )  -  qu  =  0 

with  the  boundary  conditions 

(3.hk)  .B^(u)  =  BgCu)  =  0. 

Consequently,  w^Cx)  is  an  eigenfunction  of  the  operator  L  with  the  boundary 
conditions  (3.4^*),  and  the  corresponding  eigenvalue  is  X  =  0, 

The  preceding  discussion  may  be  summarized  in  the  following 
Theorem  VI.  Let  w. (x)  be  the  solution  of 

Lw  =  -(pw  ;  +  qw  =  0 

satisfying,  the  unmixed  boundary  condition  B  (w)  =  0.  and  let  v^Cx)  be  the 
solution  of  the  same  differential  eqixation  satisfying  the  unmixed  boundary 
condition  E  (w)  =  0;  then  the  Oreen's  fujiction  of  L  with  the  boundary  condi- 
tions B  and  B^  is  given  by  the  following  formula; 

w^(x)w2(t)H(t-x)  +  W2(x)w^(t)H(x-t) 
e^^'*^  =   J(w2.w^) 

where  J(w^.w  )  is  the  conjunct  of  w^  and  w  evaluated  at  any  point. 
Problems ;  P 

3.15  Find  the  Green's  function  for  L  =-  — r  in  the  manifold  defined  by  the 

boxmdary  conditions  u  (o)  =  0,  u(l)  =0.    _ 

3.16  Find  the  Green's  function  for  L  =  -  — r  -  k  ,  k  constant,  in  the  manifold 

dx*^ 
defined  by  the  boundary  conditions  u(0)  =  0,  u(l)  =  0.  For  what  values  of  k 

doee  the  formula  break  down? 
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3.17     Find  the  Green's  fmiction  f or  L  =  -  -^  (x  t^)  , 


I 


n  the  iDanifold  defined  hy 


the  "boTindary  conditions  u(0)  =  0,  u  (l)  =  0, 

3.18  Use  the  Green's  ftmction  to  find  e  function  nix)   such  that  u  +  k  u  =  f (x) 
and  such  that  u  (o)  =  u  (l)  =  0. 

3.19  Show  that  fornnils   (3»50)  satisfies  the  conditions  of  Theorem  V, 

The  Case  where  the  Homogeneous  Equation  has  a  Non-Trivial  Solution 

Suppose  there  exists  a  non-zero  solution  w(i)  of  the  homogeneous  equation 


sxjch  that 


-  Lw  =  (pw  )  -  qw  =  0 

XX 


\iv)     =  BgCw)  =  0. 


Consider  now  the  problem  of  finding  u(x)  such  that 


with  the  conditions 


(pu^)^  -  qii  =  f  (x) 


B^(n)  =  BgCu)  =  0. 


It  can  he  proved  that  if  L  is  self-ad.ioint  end  hns  a  closed  range  and  if  the 
homogeneous  equation 

La  =  0 

has  a  single  non- trivial  solution  w,  then  the  non-homogeneous  equation 

La  =  f 

has  a  solution  if  end  only  if 

(f.w)  =  0. 

We  shall  show  that  a  similar  theorem  is  valid  for  differential  operators. 
Theorem  711.  Let  w(x)  he  the  solution,  tmigue  except  for  a  cons tent  factor, 

of  the  homogeneous  differential  eq-gation 

/   '>  • 

Iw  =  -Vpw  ;  +  qw  =  0 

such  that  B  (w)  =  B  (w)  =  0;  then  the  non-homogeneous  equation 

(3.51)  Lu  =  f 

has  a  solution  satisfying  the  boundary  conditions  B  (u)  =  Bp(u)  =  0  if  end 
only  if  ^ 

(3.52)  f      f (x)w(x)dx  +  0. 

We  shall  here  assume  that  the  hoxindary  conditions  ere  unmixed.  The  theorem 
is  valid  for  arbitrary  boundary  conditions. 

If  a  function  u  exists  satisfying  (3.5l)  and  the  boxindary  conditions,  then 

since  L  is  self-adjoint,  we  have 

0  =   (u,Lw)  -  (Lu,w)  =  -(f.w) 
and  this  shows  that  (3.52)  is  satisfied.  Conversely,  if  we  suppose  that  (3.52) 
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is  satisfied,  then  we  may  solve  (3.51).  Let  v(x)  be  a  solution  of 

and  let  t(x)  "be   independent  of  w(x).  Then^'by  using  the  Green's  function  tech- 
nique as  illustrated  in  Example  One,  we  find  that  a  solution  of  (3.51)  which 
satisfies  the  conditions  u(0)  =  u'(0)  =  0,  instead  of  the  conditions 
B^(u)  =  B2(u)  =  0.  18 

(3.53)  u(x)  =  ^  /""   fvdT  -  -^  A  fV  dT  . 

p[w,vj  -^  pb'."»'J  ^0 

Note  that,  as  shown  pre-riously,  p[w,v]  is  constant  for  all  values  of  x. 

We  can  now  show  that  the  function  defined  by  (3.53)  will  "be  the  solution 
of  the  original  prohlem.  It  is  clear  that  (3.53)  satisfies  the  differential 
equation  (3,51).  We  now  Investigate  the  houndary  conditions.  Since  the  bound- 
ary conditions  are  linear,  homogeneous,  and  unmixed,  we  may  write  them  as 
follows: 

B^(u)  =  a^u(O)  +  P^u«(0) 

B^M     =     agud)  +  Pg^'d). 

where  a.  ,P- ,a_,   P^,  are  given  constants.     Since  u(x)  in  (3«53)  was  constructed 
80  that  -0(0)  =  u'(0)  =  0,   it  is  obvious  that  B.(u)  =  0.     From  (3.53)  we  have 

(3.51*)  u(l)     =     -^"    /      fv  d'^   . 

p[w,vj   (/O 


because  from  (3.50), 

A 

fw  d  T    =     0. 


/ 


If  we  differentiate    (3.53)  usin^j  condition  (3.52)  end  put  x  equal  to  imity,  wp 
obtain 

(3.55)  u»(l)    =    ^^    /^  fv  dT. 

pLw,v]    c/o 

Fonmilea   (3.5^+)  and   (3.55)  together  show  that 

B„(w)      fl 

B  (u)   =  -4—  /     fv  d  r  =   0 

because  of  the  aasujiiption  on  w«     We  have  therefore  proved  Theorem  VII  and  we 
have  also  shown  that  if  (3.52)   is  satisfied  then  the  solution  of   (3.51)  which 
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satisfies  the  conditions 

Ti(0)  =  u«(0)  «  0 
will  also  satisfy  the  "boundary  conditons 

B^U)  =  B^Cu)  =  0. 

We  could  not  use  the  Green's  ftinction  technique  in  this  case  'becaxise  the 
equation 

(3.56)  (pu   )x  -  q^l  =  -   6(x-t) 

with  the  "boundary  conditions  B   (u)  =  B   (u)  =  0  has  no  solution,   since 

y 6(x-t)w(x)dx  /  Oj 

that  is,  (3.53)  is  not  satisfied.  If  the  right  side  of  (3.56)  is  modified  so 
thst  it  Is  orthogonal  to  w(x),  vre  o"btaln  a  pro"blem  which  can  "be  solved.  This 
is  done  "by  writing 

(-5.57)        ^P^x^x  -  '^''^  =  -   ^(^-*)  ■*•  ^^Mt 
where  V  is  found  from  the  equation 
w(t)  -  vy  w^dx  •  0. 

(■^.57)  can  "be  solved  in  the  same  •.^ray  as  (3. 51)  was.  The  details  will  be  left 

to  the  reader. 

Problem: 

3.20.   In  each  of  ttie  following  cases  find  the  appropriate  orthogonality  relation 

that  f(x)  must  satisfy  in  order  that  the  equation  may  have  a  solution,  and  then 

find  the  solution: 


(a)  u"=  f(x).  u'(0)  =  u'(l)  =  0; 

(1))  u"+  TT^u  =  f(x),  u(0)  =  u(l)  =  0; 

(c)  (ru')'=  f(x).  xu'(x)|^Q  =  0,  u'(l  )  =  Oi 

(d)  ^i*   =  fU).  u(0)  =  u(l),  u'(0)  =  u'(l). 
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Green 'g  ''^unction  for  (General  Boimdary  Conditions 

In  the  penersl  case  where  the  'boundary  conditions  ere  mixed,   the  Green's 

f-unction  may  still  "be  found  in  a  straightforward  way.     Suppose  the  equation 

ie 

(:>.'?8)  (pr   )     -   q^  =  -  8(x-t) 

X    X 

with  the  "boundary  conditions 
(3.59)        \(p)   =  BgCf-)  =  0. 
I«t  v(x)  and  w(x)  "be  any  two  lineprly  independent  solutions  of 


(pu^)^  -  qu  =  0. 


Then  write 


(-^.60)    gU.t)  =   a,v(.)  ■.  P,w(x)  +  v(x)w(t)H(t-x)-Kw(x)v(t)H(x-t)^ 

^        ^  Cw.^  p(t) 

where  o-.  and  P,  are  constants  which  will  be  determined  so  that  g(x,t)  satisfies 
the  houndery  conditions  (3.59) •  -^he  "bracket  expression  [w,vj  is  the  Wronsklan 
of  w  and  v  at  the  value  t.  Note  that  it  cannot  he  zero  for  any  value  of  t 
"because  if  it  were  zero,  then  hy  the  Corollp-ry  to  Theorem  IV  and  hy  Theorem  II. 
the  functions  v(x)  and  w(x)  would  not  he  linearly  independent. 

We  show  first  thet  g(x,t)  satisfies  the  appropriate  conditions  at  x  =  t, 
namely  thfit  ^(x.t)  is  continuous  and  its  derivative  has  a  jump  of  magnitude 
-l/p(t).   The  continuity  i<3  ohvious.  From  the  dpflnitlon  of  the  Wronsklan  it 
follows  that  the  Jump  in  the  derivative  at  x  =  t  is 

w'(t)v(t)  -  v'(t)w(t)  _     1 

Ew.v:]p(t)  ="?(Ty  • 

U'ow  conslde-^  the  houndnry  ccnditions.  Since  they  are  liner^r  and  homo- 
geneous ,  we  have 


(3.6])        B^(g)  =  a^B^(v)  +  P^Bj^Cw)  +  B^(r) 


BgCe)  =  ^^^^(v)  +  \^2^^)   +  B2(r), 
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where  we  have  used  r(x)  to  represent  the  function 
v(x)v(t)H(t-x)  +  w(x)v(t)H(x-t) 

&,£]p(t) 

The  "boTindary  conditions  (3.59)  will  consequently  give  us  two  linear  equations 

to  determine  the  values  of  cl,  and  P, .  These  equations  will  have  a  solution  if  the 


determinant 
(3.62) 


B^(v)    B^(w) 


B  (v)    B  (w) 


does  not  vanish. 

If  (3.62)  vanishes,  then  either  B^(v)  =  B^{y)  =  0,  which  implies  that  v 
is  an  eif:enfianction  of  the  operator  corresponding  to  the  eigenvalue  eero,  or 
there  exists  a  constant  c  such  that 

Bj^(w)  +  c  B^(v)  =  0 

B2(w)  +  c  ^^W)   =  0. 

These  equp.tlons  Innily  that  B^(w  +  cv)  =  B   (w  +  cv)  «  0;  hence  w  +  cv  is  an 
eigenfunctlon  of  the  operator  and  corresponds   to  the  eigenvalue  zero. 


V/e  see   then  that  if  there   is  no  eigenfunctlon  corresponding  to  the  eigen- 
value zero,   the  Green's  function  exists  and  is   given  by  (3.60)  where  a  and  p 
are  found  by  solving  ( 3 . $9 ) . 

Green's  Fiinction  for  an  n-th  Order  Equation  or  a  System 

The  method  of  the  preceding  section  can  be  easily  generalized  in  order  to 
obtain  the  Green's  function  for  an  n-th  order  differential  equation  or  for  a 
system  of  n  first-order  differential  equations. 

The   Green's  function  for  an  n-th  order  differential  operator  is  defined 
to  be  the  symbolic  'function  g(x,t)  which  satisfies  the  equation 

(3.63)  Lg  =  p     2_£     +  p     2. i  +...  +  P         5^+pg  =  -5(x-t) 

^■^     -''  ^       ^o   ^  n  1  J  n-1  ^n-1  dx       ^n^ 

dx  dx 

and  the  boundary  conditions 


\{%)   =  \k%)   =  ...  -  Bjg)  =  0. 
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Since  the  8-function  is  zero  except  when  x  =  t,  equation  (3«63)  implies  that 
g  satisfies  the  homogeneous  equation 

Lg  «  0 
for  all  values  of  x  other  than  x  ■  t.     If  we  integrate  (3.36)  with  respect  to 
X  from  0  to  X  we  get 


^n-1 

o    ,  n-1 
dx 


n  1 

_  /  p.lJLs  dx+/ 

V     ^o    ,  n-1  » 

o         0  dx  o 


^n-1 

Pl-lrf*  — *PnS 
dx 


dx  -  -H(x-t). 


If  we  assume  that  the  functions  Pi^(k  =  0,  l,,n)  have  continuous  derivatives  of 


n-1 


order  n-k  then  this  result  implies  that  p  ?  has  a  jmnp  of  magnitude 

(-1)  at  X  =  t»  Further  integrations  show  that  the  functions  p 


d   g 
3    '  ■  z   ,  ,,,  p  g  are  continuous  at  x  ■  t, 
o  ,  n— 3       o 

dx 


dx^ 


Assume  that  p   (x)  is  continuous  and  non-zero  in  the  interval  (0,1) j   then 
we  conclude  that  if  g  satisfies  (3.63)  it  will  also  satisfy  the  following  con- 
ditions: g(x,t),  together  with  its  first  n-2  derivatives,  will  be  continuous  at 
X  =  t,  while  the  (n-l)th  deriva.tive  will  have  a  jump  of •  magnitude  -l/p  (t) 
there.     This  conclusion  caji  be  easily  verified.     Put 

_  _vn-l 
(3.61;) 


^^'''^^  "  "     (i-in  p  (t)     "^"""^^  *  k(x,t),- 


then  g  will  satisfy  the  above  conditions  if  k(x,t)  and  its  first  (n-l)  deriva- 
tives are  continuous  at  x  =  t.  Now  if  (3»6U)  is  substituted  in  (3.63)  we 
find  that 

(3.65)  Lk  -  r(x,t), 

where  r(x,t)   is  a  piecewise  continuous  function  of  x.     Just  as  in  Appendix  II, 
we  may  prove  that  equation  (3.65)  has  a  continuous  solution  such  that  the 
first  n-1  derivatives  of  this  solution  is  lite  wise  continuous. 

Now  since  the  solution  of  (3.63)  which  satisfies  the  boundary  conditions 
is  unique,  g(x,t)  must  be  of  the  form  (3.6U)  and  therefore  must  obey  the  con- 
tinuity conditions  stated  above. 
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The  definition  of  the  Green's  function  for  a  system  is  more  complicated, 
and  will  be  considered  In  Appendix  I. 

Green's  yTinction  and  the  Boundary  Value  Problem 

So  fer,  we  have  discussed  the  solution  of  non-homof eneous  differential 
equations  with  homogeneous  boundsry  conditions.  We  now  wish  to  consider  the 
case  of  non-homog-eneous  boundary  conditions,  both  for  homogeneous  and  non- 
homogeneous  differential  eqtiations.  We  shall  show  that  the  Green's  function 
technique  can  be  used  to  solve  these  more  conplicated  probleoss  also, 

?or  example,  suppose  thst  we  wish  to  find  a  solution  u  of  the  differential 
equation 

(-'.66)  La  =  if  +  a(x)u'  =  f 

such  that  u(0)  =  b^,  u(l)  =  b  .  Let  h(x,t)  be  the  Green's  function  of  the  adjoint 
homogeneous  eqriation;  then 

(3.67)       ^*^  =  H  -  ^  ^^^^  =  ^^^-*^ 

dz 

and  h(x,t)  satisfies  the  homogeneous  adjoint  boundary  conditions 

h(O.t)  =  h(l,t)  =  0. 

Mi:!.ltiply  equation  (3.66)  by  h,  equation  (3.6?)  by  u,   subtract,  and  Integrate 

over  (0,l),  We  get 

11  1 

/u(x)6(x-t)dx  -yh(x,t)f(x)dx  =  <u,L*h)  -  (h.Lu)  =  (uh'-hu') 

0  0 

=  u(l)h'(l,t)  -  u(0)  h'(0,t) 

or,  finally, 

1 
(^.68)    u(t)  = /h(x,t)f(x)dx  +  b-h'd.t)  -  b  h'(0,t). 

•^  1  'O 

o 
The  primes  denote  differentiation  with  respect  to  x. 

The  same  method  can  be  used  In  more  general  cases.   "Before  proving  this  fpct, 
we  must  investigate  a  little  more  closely  the  boundary  conditions  of  the  adjoint 
problem.  'From  eoustions  (3.9)  and  (3.10)  we  have 

1 

(3.69)       ■  (r.Lti)  -  (u,L*v)  =   J  (v,u) 

o 

Suppose  that  the  boxmdary  conditions  associsted  with  the  operator  L  are 

B^(u)  =  b^,  B2(u)  =  bg.-",  B^(u)  =  b^. 

Since  the  boundary  conditions  depend  linearly  on  the  2n  values  of  u  and  the  first 
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(n-l)  derivatives  of  u  at  x  =  0  and  x  =  1,  we  may  solve  for  any  n  of  these 

najn"bcr8,  say  u(o),  u  (o),,..,  u^  '(O),  u(l),  u  (l),...,  u^^"  ~     ,  in  terms  of 

!}.»•.•»  Tj  *     When  these  expressions  are  substituted  into  J(v,u)^   ,  we  get 
in  o 


(3.70) 


J(v.u) 


=    h.C-(v)+..»  t  C   (v)  +  other  terms. 


A  simple  calculation  shows  that  the  expressions  C  (v),...,  C  Cv;  are  linear 
homogeneous  expressions  containing  the  values  of  v  and  its  first  (n-l) 
derivatives  at  x  =  0  and  x  =  1,  The  "other  terms"  in  (3»70)  are  linear 
homogeneous  expressions  involving  the  n  ntanhers,  u   (o),.,,,  u     (o), 
u^'^'^'^^d ),...,  ti^°^^\l).  We  shall  denote  these  n  numhers  V 
C.(u),,,.,  C  (u)  respectively;  then  the  "other  terms"  in  (3»70)  may  he 
written  as 

-  C  (u)  B  (v) C  (u)  b"""  (v) 

1     1  n     n    > 

where  B^(v) , • • • ,B^(v)  are  linear  homogeneous   expressions    involving  the  values   of 
V  and  its   first    (n-l)   d-^rivatlves  at  x  =   0  and  x  =   1.     Finally,    Instead  of    (3.70) 
we  have 

(•^.71)  J(u,v)   ^  =  B   (u)C*(v)+   •..  +  B   (u)C*(v)  -   C,  (u)B*(v)_   ..._C   (u)B*(v); 

o         ±  1.  nn  IX  nn 

consequently,  the  homopeneons  adjoint  "bo-'incJr'.ry  conditions  are 

B  (v)  =  ...  =  B*(v^  =  0. 
•!•  n 

As  an  illustration,  suppose  that  Lu  =  u     nnd  that 

B^(n)  =  u(0)+u(l),  B2(u)  =  u'(o).  B^(u)  =  u"(0)  +  u'(l).  B^(u)  =  u"(l). 

The  operator  1  Is  formnlly  self-adjoint,  and  "by  evaluating    (v.Lu)  -    (u,Lv)  we  find 

that 


J(v,u) 


Eliminating  u(0),  u'(o),  u'(l)  and  u"(l)  from  this  formula,  we  g^t 
J(v,u)  I   =  v(l)u"«(l)-v'(l)Bjj^+Tr"(l)(B^-u"(0))-v'«(l)u(l) 
-  v(0)u"'(0)+v'(0)u"(0)-v"(0)B2+v"«(0)(B^-u(l)) 
=  B^v"'(0)-B2v"(0)  +  B^v"(l)  -  B^v'(l) 

+  u«"(l)v(l)+u"(o)(v'(0)-v"(l))-u"'(o)v(o)-u(l)(v"«(l)+v"'(0)); 
consequently, 

B^(v)  =  -v(l),  B*(t)  =  -v'(0)+v"(l).  B*(v)  =  v(0),  B*(v)  =  v '"  (l)+v"«  (o). 
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Let  us  now  return  to  the  'bound.ary  vedue  problem.  Suppose  that  the  differential 
equation  to  "be  solved  is 

(3.72)  Lu  =  f 

with  the  non-homogeneouB  boundary  conditions 

B   (u)     =     h    ,...,  B   (u)     =     b   . 
1  in  n 

Let  h(x,t)  be  the  Green's  frmction  of  the  adjoint  homogeneous  equation:  thus 

(3.73)  ^*h.     =  6(x-t). 

and  h(x,t)  satisfies  the  homogeneous  adjoint  boundary  conditions 

B*(h)  =  ...  =  B*(h)  =  0. 
1  n 

Multiply  equation  (3.72)  by  h,  equation  (3.73)  by  f ,  subtract  ,  and  integrate 
over  (0,l).  Then,  using  (3.7l)  we  get 

f   u(i)6(x-t)dx  -  /*  h(i.t)f  (x)dx  =  (u,L*h)  -  (h.Lu)  =  -  J(h,u)  | 

=  B,(u)C*(h)+  ...  +B  (u)C*(h)-C.  (u)B*(h)-  ...  -C  (u)B*(h). 
11         nnli         nn 

If  we  make  use  of  the  boundary  conditions  on  u  and  h,  this  expression  reduces 
to  the  following: 

(3.7^)  u(t)  =  /*  h(x,t)f(x)dx+  b  C  (h)  +  ...  +  b^C^(h). 

This  proves 

Theorem  7 I II,  The  solution  of  Lu  =  f  which  satisfies  the  boundary  conditicns 

B  (u)  =  b- \^^^  -   \  ^^   given  by  (3.7^)  where  h(x.t)  is  the  Green's 

function  of  the  ad.^oint  homoisceneous  equation  (3.73). 
4<        ■I' 
Note  that  C  (h),...,C  (h),  which  appeer  in  (3-7^),  are  actually  functions 

of  t.  Note  also  that  the  method  indicated  above  works  just  as  well  for  a 

system  of  first-order  differential  equations  as  for  an  n-th-order  differential 

equation. 

Problems: 

3.23  Solve  u  +  u  =  f (x)  with  the  conditions  u  (O)  =  0,  u(o)  +  u(l)  =  2. 

3.2i+  Solve  u     =  f  with  the  conditions  B  (u)  =  1,  B-(u)  =  0,  B  (u)  =  0 

and  B^(u)  =  1,  where  B  jB^.B^  end  Br  are  the  expressions  used  in  the  text 

in  connection  with  this  operator. 
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3.25  Solve   the  system 

with  the  "boTindary  conditions  u_  (o)  =  1  and  u^d)  =  2, 

(Hint,  Refer  to  Appendix  I  for  the  Green's  matrix  of  a  solution). 

Crreen's  Function  of  the  Ad.joint  Equation 

Suppose  that  in  Theorem  VIII  we  consider  the  case  of  homogeneous  "boundary 
conditions,   that  is,  b.   =   ...  ^  h     =0;    then  the  solution  of  liu  =  f  as  given 
hy  (3.7^^)  is  ^ 

u(t)  =  r     h(x,t)f(x)di. 

However,  we  know  that  the  solution  can  also  he  written  as 

.1 
u(x)  =  r       g(x,t)f(t)dt, 

0 

where  g(x,t)  is  the  Green's  function  of  the  original  operator  L.  This  indicates 

that  there  is  a  connection  "between  g(x,t)  and  h(x,t).  We  prove 

Theorem  IX,  Let  g(x.t)  be  the  Green's  function  for  an  operator  L  on  a  manifold 

defined  by  certain  boundary  conditions,  and  let  h(x.t)  be  the  Green's  function 

for  the  ad.joint  operator  L  on  the  manifold  defined  by  the  ad.joint  boundary 

conditions;  then 

g(x,t)  =  h(t,x). 

In  particular,  if  L  is  self-ttd.joint  then  g(x.t)  is  a  symmetric  function  of  x 

and  t. 

To  prove  the  theoreia  we  start  with  the  defining  relations  for  the  Green's 

functions,  namely 

l€(x,t)  =  6(x,-t), 

(3.75)  , 

Lh(x,Y)  =  6(x-r). 

Prom  the  properties  of  the  adjoint  operator,  we  get 

yh(x,Y)l€U,t)dx  =yL*h(x.r)g(x,t)dx. 

If  we  use  (3.75),  this  equation  becomes 

yh(i,T)6(x-t)dx  =y6(x-T)g(x,t)dx, 

h(t,r)  =  g(r,t)j 

this  last  result,  with  a  slight  change  of  notation,  is  the  conclusion  of 

Theorem  IX. 

Problem; 

3.26  Verify  Theorem  IX  if  1  =  -i- +  x  -L,  with  the  conditions  u(0)  =u(l)  and 

'  ^  V    '  ,  V  J  ^     dx 

u  (0)  +  u  (l)  =0.  °^ 
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Discontinuity  Conditions 

In  many  practical  applications  it  is  important  to  consider  linear  differ- 
ential eq-uations  in  which  the  coefficients  may  "be  discontinuous.  This  situa- 
tion occurs  in  cases  where  the  properties  of  the  medium  vary  discontlnuoasly, 
for  ejca-nple  the  propagation  of  sound  waves  through  media  of  different  densi- 
ties or  the  tiansmission  of  heat  through  materials  of  different  thermal  conduc- 
tivity. The  differential  equation  describes  the  "behavior  of  the  solution  in 
each  medium  separately,  "but  it  is  still  necessary  to  determine  how  the  solution 
"behpves  pcross  the  Interface  of  the  two  media.  Usually  the  behavior  of  the 
solution  across  the  inte^-face,  from  which  we  derive  the  so-called  discontinuity 
cnndltlonR.  is  determined  "by  physical  considerations;  for  example,  in  sound 
propapption  the  pressure  and  velocity  of  the  sound  wave  mast  be  continuous 
across  the  interface,  otherwise  there  would  be  an  infinite  acceleration  there. 
Slmllnrly,  in  electromppnetic  theory  we  can  show  that  because  of  Ga-uss'  theorem 
the  tfiTventlel  components  of  the  electric  field  must  be  continuous  across  the 
interface.  Clerrly,  these  discontinuity  conditions  are  as  much  a  part  of  the 
physical  description  of  the  problem  as  are  the  differential  equations  themselves 
and  consequently  the  question  of  discontinuity  conditions  is  not  a  mathematical 
question.   However,  it  is  useful  to  investigate  what  kind  of  discontinuity 
conditions  sre  mathematically  appropriate  for  given  differential  equations. 

Consider  the  self-adjoint  differential  equation 

(3. ■'6)         -lAi  =  (pu')'  -  qu  =  f(x) 

where  p(i),  q(x),  and  f (x)  are  continuous  except  for  possible  jumps  at  the 
points  X  =  i^,  X2.*",Xj..  Theorem  II-I  tells  us  that  despite  the  jumps  there 
exist  solutions  of  (3. '^6)  such  that  the  functions  u(x)  and  p(x)u'(y)  are  con- 
tinuous throughout  the  entire  interval; consequently  we  may  write  the  discontin- 
uity conditions  for  this  equation  as  follows: 

u(x  +0)  =  u(x.-O) 

T)(xj+0)u'(xj+0)  =  p(xj-0)u'(x^-0),   i  =  1,2, ...k. 

T'fote ,  that  if  p(i)  is  continuous  st  x  =  x. ,  the  last  condition  states  that 
n'(x)  must  be  continuous  at  x  =  x . . 

Instead  of  using  the  reaiilts  of  a  theorem  such  as  Theorem  II-l  about  the 
existence  of  a  sol-'Ttlon  of  the  differential  equation,  we  may  obtain  the  dis- 
continuity conditions  by  the  following  simple  argument.   Integrate  equation 
(3.^6)  across  a  discontinuity,  say  from  x  =  x^-c  to  y  =  xj+g,  where  e  js  a 
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small  Tiosltlve  quantity.     We  have 


pu' 


+  f )dx. 


Since  q,  u,  and  f  are  integra"ble  ftinctions,  the  right-hand  side  of  this 
equation  apnronches  7ero  as  e  approaches  zero,  and  therefore  pu'  is  contin- 
uous across  x..     Now  write  (3.76)  in  an  intef:ral  fonS: 


"■  =  "^m^  ^  ;b/<-->^^  ' 


This  is  valid  if  p(x)  ^   0.   If  we  integrate  again  across  x,  we  get 

x.+c     X.+€  x.+c 

u     =/     2iO)ulWL^^J  dx   /(^,,j,^  ^ 


The  integrals  again  approach  zero  as  e  approaches  zero  and  therefore  u(x)  also 
is  continuous  at  x  =  x. . 

The  same  concltieinns  about  the  continuity  of  u(x)  and  pu'  may  he  obtained 

different  axi 
then  we  can  vrrite 


"by  a  different  argument.  Suppose  u(x)  has  a  Jtunp  of  magnlttuie  a  at  i  =  x  ; 


u(x)  =  aH(x-Xj^)  +  v(x), 

where  v(x)  is  continuous  at  x  =  x. •  We  would  now  find  that  Im.  is  not  a  piece- 
wise  continuous  function  "but  instead  contains  the  derivative  of  a  delta  func- 
tion.  Similarly,  if  pu  has  a  Jnmp  discontinuity,  Lu  contains  a  delta  function. 

Since  Lu  is  given  as  a  piecewise  continuous  function,  it  follows  that  both  u  and 

I 
pu  must  "be  continuous.  We  state  our  conclusions  in  the  following 

t 
Rule;  The  discontinuity  conditions  for  the  operator  Lu  =  -(pu  )  ■♦•  qu  are  that 

I 
both  u  and  pu  "be  continuous. 

This  rule  has  an  important  consequence  which  is  based  on  the  concept  of 

the  impedance  of  a  solution.  We  define  the  impedance  of  a  solution  u(x)  of 

/   '\  Pu' 

the  equ8.tion  La.  =  -(pu  )  -  qu  as  the  ratio  •  Then  the  above  rule  implies 

that  the  impedance  is  continuous  across  a  discontinuity, 

Aa  an  application  of  these  rules  we  discuss  the  propagation  of  a  plane 
wave  through  two  media  of  different  properties.  The  mathematical  problem 
that  corresponds  to  this  physical  question  is  the  following: 
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Consider  the  equation 

'    2 
u  +  k  u  =  0,  -00  <  X  <  00 

where  k  has  one  constant  value,  k. ,  for  x  <  x  and  has  another  constant 
valtie,  k^,   for.x  >i  .  Find  a  solution  u  of  this  equation  such  that,  for 

X  very  large  nef^etively,  11=6'^  . 

ikiX 
For  x<  Xj,,e  1  is  a  solution  of  the  differential  equation.  For  x  >x  , 

^ik2(x*Xo)  ^  ggik2(x-x^;  ^^   ^^  ^^^^  ^^  ^^^  general  solution  of  the  equation. 

Since  the  properties  of  the  medium  change  discontinuously  across  x  =  x  ,  then 

I 
by  the  nile  stated  above,  both  u  and  u  (note  p  is  identically  unity)  must  be 

continuous  at  x  =  x  *  Diese  conditions  give  the  following  equation  for  A  and  B: 

0 

e^(ik^k^)  =  A  +  B 


ikj^exp(-ik^x^)  =  ikgCA  -  B), 


By  division  we  find  that 


and  then 


-.i 

-  B 
+  B 

=     h 

A 

^2 

*^ 

B 

^2 

-^1 

The  fact  that  this  result  could  have  been  obtained  directly,  by  using  only  the 
Impedance  rule,  illustrates  the  advantages  of  this  rule.  It  enables  us  to 
obtain  the  ratio  A/B  inmedlately  without  solving  a  set  of  linear  equations. 
In  more  conrpllcated  cases  where  there  is  lore  than  one  discontinuity  present, 
such  a  use  of  the  impedance  rule  may  produce  a  great  simplif icstion. 

Of  course,  to  find  the  values  of  A  and  B  we  can  solve  the  original 
equations.  Die  solution  is 

^=  -2kJ"  «^(iVo^'  ^=  "a^«^^-^Vo^- 

The  physical  interpretation  of  the  solution  is  that  a  plane  wave  Beip[;-ik2(x-x  1 
comes  In  from  the  right;  one  part,  namely  A  expLlk_(x-x  )"]  ,   is  reflected  at  the 
interface;  the  other  part,  exp(-ik-x),  is  transmitted  across  the  interface.  In 
a  later  chapter  we  shall  discuss  these  questions  more  thOTOughly. 

It  should  be  noted  that  caution  is  needeed  in  the  application  of  the  above 
rule  for  discontinuity  conditions.  After  all,  the  discontinuity 
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conditions  sre  psrt  of   the  physical  problem  and  cannot  be  obtained  by  mathe- 
maticpl  considerations  alone.  However,  if  the  mathematical  eqiiations  are 
properly  formulated  the  rule  will  give  the  correct  discontinul-ty  conditions. 

To  understand  what  is  meant  by  a  proper  formulation^  consider  the  case  of 
one-dimensional dectromasnetic  wave  propagation.  In  general,  two  kinds  of 
propp^ation  are  possible?  they  are  called  transverse  electric  (TE)  and  trans- 
verse magnetic  (TM)  propagation.   In  both  TE  and  TM  propagation  we  must  solve 

the  equation 

2 

u"  +  It)  ^l€  U  =  0, 

where  oi  is  the  frequency  of  the  ware,  \j.   the  magnetic  permeability,  and  e  the 
dielectric  constant  of  the  medium.  Suppose  we  are  interested  in  the  trans- 
mission of  a  wave  across  the  interface  of  two  media  for  which  e  has  a  constant 
value,  but  different  in  each  medium,  while  the  constant  value  y.   is  the  same 
for  both  media.  IPor  IB  propagation  the  discontinuity  conditions  are  that 

u  and  u  be  continuous  across  the  interface,  while  for  TH  propagation  the  dis- 

I 
continuity  conditions  are  that  u  and  en  be  continaotis.  Zhia  later  case  seems 

to  contradict  our  rule.  However,  if  we  scrutinize  the  physical  situation  we 

observe  that  for  the  TE  case  the  eqtiatlon 

2 
u"  +  0)  tieu  =  0 

is  correct  even  when  c  varies  continuoxosly,  but  for  the  TM  case  the  correct 

equation  if  e  varies  continuously  is 

1  ? 

-  (cu')'  +  u)  cM-u  =  0. 

e 

If  we  apply  the  nale  to  this  latter  case  we  get  the  correct  discontinuity 
conditions,  namely  u  and  cu'  continuous. 

"Prom  this  illustration  we  deduce  the  requirement  that  the  mathematical 
eauPtlong  should  be  so  formulated  that  they  pre  valid  for  continuous  changes 
of  the  pprameters  specified  in  the  problem..   If  they  are  so  formulated,  then 
the  r^ile  of  this  section  will  give  the  physically  correct  discontinuity  conditions. 

We  shall  give  one  final  e»mple.  Consider  the  equntion 

u"  +  a.   6(t-Xjj)  u  =  0, 

Wh».t  is  to  be  understood  by  this  equation?   In  order  that  the  term  u(x)  times 
the  delta  function  malce  sense,  u(x)  must  be  continuous  at  x  =  x  ,   Integrate 

0 

the  equation  from  x  -c  to  x  +€.  We  get 

0         0 


u' 


0 
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therefore  the   jump  in  u'(x)   at  x  =  x    must  equal    -att(x  ).     For  values  o£  x  ^  x  , 
the  value  of  the  delta  function  is  zero,   and  the  equation  reduces  to 

u     =  0    . 

We  conclude  that  the  equation  is  to  be  interpreted  as  follows: 

The  function  u  is  a  solution  of  u"=  0  for  all  values  of  x  different  from  x  . 

o 

At  X  =  x  ,  the  function  u  is  continuous,  but  its  derivative  has  a  jump  of 

magnitude  -eu(x  ) .  With  these  specifications  the  solution  of  the  equation 
is  completely  determined  as  soon  as  two  boundary  conditions  are  given. 

Problems  ; 

3.27.  Find  the  Green's  function  for  Lu  »  k  u  -  u"  (-oo<  x  <  oo  )  when  k  is  a 
constant  k^  for  x  <  0,   another  constant  k_  for  x  >  0  and  when  the  boundary  con- 
ditions are  that  u  vanish  at  both  +  go, 

3.28.  Find  a  solution  of  u"  +  k  u  =  0  such  that  u  «»  e       -^     for  large  negative 
values  of  x,   given  that  k  is  a  constant  k,   for  x  <  0,   another  constant  k„  for 
0  <  X  <  a,   and  again  the  constant  k^   for  x  >  a, 

3.29.  Find  u(x)   satisfying  u"  +  a6(x-x^)  u»  0  . 

Singular  Points  of  a  Differential  Equation 

During  the  previous  discussion  of  differential  equations  we  have  always 
made  the  assumption  that  p(x)  j^  0,  because^  ss  was  pointed  out  in  an  earlier 
section,  if  p(x)  «  0  then  some  of  the  solutions  of  the  equation 

-(pu')     +qu=0 
may  become  infinite  at  the  points  where  p(x)   =  0.     However,   since  in  many 
physical  applications  p(x)   does  become  zero  at  some  point  in  the  interval, 
it  is  important  to  investigate  the  behavior  of  the   solutions  of  the  differen- 
tial equation  at  such  points. 

In  this  section  ve   shall  discuss  the  singularities  of  differential  equa- 
tions and  determine  the  nature  of  the   solutions  of  the  homogeneous  equation  at 
such  points   .     In  a  later  section  we   shall  investigate  the  behavior  of  the 
solutions  of  the   inhomogeneous  equation.       We   s'hall  find  that  just  as  in  the 
theory  of  power  series,   in  order  to  give  a  complete  treatment  of  the  nature 
of  the  singularity  we  must  consider  the  equation  and  its  solutions  not  only 
for  real  values  of  x,  but  also  for  complex  values. 
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Consider  then  the  n-th-order  differential  equation 

p^(x)u^''^(x)   +  p^(x)u^""^\x)   +...+  p^_^(x)u'(x)  +   p^(x)u  =  0, 

v;here  p    ,p,  , . . .  ,p     are   single-valued  analytic  functions  of  x  in  some  region 
0     1  n 

R  of  the  complex  plane.     Rewrite  the  equation  es  follows: 

u^""^  +  P-^/p^u^"-^^    ...  +  P^.i/PqU'    +  pyp^u  .  0    . 

Any  point  x  at  which  all  the  function  p  /p   .p^/p    , . . .,  P^/Pq  ^^  regular, 

that  is,   have  no  singularities,  is  called  an  ordinary  point  of  the  differ- 
ential equation.     It  can  be  proved  that  given  n  arbitrary  initial  conditions 
at  an  ordinary  point  x  ,   there  exists  a  function  u(x)  regular  in  the  neigh- 
borhood of  X     which  satisfies  the  differential  equation  and  the  initial  con- 
o 

ditions . 

A  point  X  at  which  at  least  one  of  the  functions  p  p^,pp/p   , ..  .,p^/p^ 

has  a  singularity  is  called  a  singular  point  of  the  differential  equation. 
For  example,   consider  the  equation 

xu"  +  u'   =  0 
which  has  x  =  0  as  a  singular  point.     The   solutions  of  this  equation  are 

u  =  a  +  p  log  X 
where  a  and  p  are  arbitrary  constants.     Vfe  notice  that  some  of  the    solutions 
of  this  differential  equation  have  a  branch-point  singularity  at  x  =  0,   the 
singular  point  of  the  equation.     The  same   situation  arises  in  the  case  of 
an  arbitrary  n-th-order  equation,     namely  some  of  the  solutions  of  the 
differential  equation  may  have  a  singularity  at  the   singular  points  of 
the  differential  equation.     All  possible  types  of  singularities  may  occur 

in  the  solutions:     poles,  branch  points,   and  essential  singularities  such 

-1/x 

as  that  of  the  function  e     '      at  x  =  0. 

We  shall  distinguish  between  singular  points  according  to  the  possible 

types  of  singrularities  the  solutions  of  the  differential  equation  may  have 

there.     A  singular  point  is  said  to  be  regular  if  all  the  solutions  of  the 

differential  equation  have   at    ioo3t    a  pole  or  a  branch  point  singularity  at 

that  point,     a  singular  point  at  which  some  solutions  have  essential  singular- 

ities  is  called  4-rre gular . 

t  For  a  more  detailed  discussion  see  \*ittaker  and  Watson , Modem  Analysis 
(Cambridge  Press)   Chnp.  VII  md  Ince,   Ordinary  Differential  Equations, 

(Dover,  19hh)   Chap.  XVI. 
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As  an  illustration  of  a  regular  singular  point,   consider  Bessel's 
differential  equation 

(3.77)  u     +iu+(l-^)u  =  0. 

X 

We   show  that  the  point  x  =  0  is  a  singular  point  of  this  equation.     Assume  that 

00 


(3.78)  u  =  x°(a^+  a^x+ag  X  +   ...)  =    JI    a^ 


n+a 

X 


> 


n=0 


where  a, a  ,a,,..,   are  constants  which  will  be  determined  to  make   (3.78) 

satisfy   (3.77).     Note  that,  we  may  assume   a  ^  0.     For  ifa=a^=...=  a,^=0 

but  a,    /  0,   then  we  may  write 

„         a+k  /      .  .  \ 

""^^      K*  Vi^-'  •••^. 

and  replacing  a+  k  by  a,   =k.  by  a  ,   a,    ^   by  a^  ,    and  so  on,  we  get  exactly  (3.78), 
We  have 

00  T  00  „ 

,       V—  c        \         n+a-1  .       x"  f    ,    \f  -i  \        n+a-2 

u'    =  >     (n+a)   ax  ,     u"  =  2_  (n+a)Cn+a-l;a  x 

n=0  "  n=0  " 

Substituting  these  expressions  in  (3.77)  we  get 

E/        N/  ^  X  n+a-2       -e—  /        \         n+a-2     ^—         n+a       2  r—         n+a-2       _ 

(n+a)(n+a-l)   ax  +  >     (n+a)a     x  +>ax       -v>ax  =0 

^^  n  ^—  n  ^—     n  ^—     n 


o  o 

or 

t   E(n*a)2-u2 J  a^x°"-2  ^  f-  ^^^n+a  ^  ^^ 

0  0 

Replace 

Z_  *  ^      ^y  Z_  a   ^3C       ; 
0  d 

then  we  have 

2_L(n+a)  -u  Ja  X     +  2_a  -x     =  0. 
0  2 

Since  this  eqiiatlon  must  hold  for  every  value  of  x,  it  follows  that  the 
coefficient  of  every  power  of  x  vanishes,  that  is, 

Z(n+a.)^-v^2   %  +  \_2   =0.  n  =  2,3,--- 


Dl+a)^-^^:!  a^  =  0 


(a^-t>^)  a  =  0. 

0 
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Since  a     ?^    0,  the  last  equation  implies  that  a  ■  1  v.     Assiune  a  =  v;  then 

the  next  to  last  equation  implies  a,    =  0  and  the  other  equations  may  be 
written  as  follows: 

\  '  -     n(n+2v)     '  ^-^2,3    ... 

This  recurrence  relation  determines  all  the  other  coefficients  in  (3.78). 

We  have 

-  a  -ap  a 

^2  "      ^.2(l+v)  '  \  '    U.2(2+v)  "   ti-2.2.2(2+v)(l+v)  '  ••* 
a^=aw=a_a...=0. 
Consequently,  the  infinite  series 

""l  =  "^  ^^  "  2.2U+V)  ■"  U.2.2.2(2+v)(l+v)  -  •  •  •  ) 
will  satisfy  (3.77)  formally.  By  the  use  of  the  ratio  test  it  can  be  shown 
that  this  series  converges  for  all  values  of  x  and  hence  it  is  actually  a 
solution  of  (3.77). 

Similarly,  if  we  assume  a  =   -v,  we  find  that  another  solution  of  (3.77) 
is  given  by  the  infinite  series 

^2  ^h 

^2'^       ^-^  '   2.291-v)  *  U.2.2.2(2-v)(l-v)   -•••)• 
Since  these  two  solutions  are  linearly  independent,  the  general  solution  of 
(3.77)  is 

u  =  avu   +  pUp 

except  when  v=Oorv=m,m  being  an  integer.  When  v  =  0,  the  infinite 
series  Up  is  exactly  the  same  as  u..  , while  when  v  =  m  the  expressions  for 

a  ,,  a  „,....  in  the  series  for  u^  become  infinite. 
m+1  m+d  d 

It  is  easy  to  handle  the  case  v  =  0.  The  equation  now  is 

(3.79)  u"  +  1  u'  +  u  »  0. 

Since   the  last  term  has  no  derivatives  and  is  regular  at  x  =  0,  we  expect 

that  the  solutions  of  (3.79)  will  behave  near  x  =  0  like  the   solutions  of 

the  equation 

"       1 
u     +  i  u'   =  0. 
x 

The  general  solution  of  this  equation  is  u  =  a  +  p  log  x;  consequently  we 
expect  that  one  solution  of  (3.79)  will  behave  like  a  constant  at  x  =  0  and 
another  solution  will  behave  like  log  x.  The  series  u,  that  we  found  pre- 
viously behaves  like  a  constant.  Now  we  must  find  a  solution  of  (3.79) 
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wliich  behaves  like  log  x.     Assume 

2  2 

Up  =  log  x(a  +  a,+  apX    +.,,)  +  (b  +  b-x  +  bpX  +   ...)    , 

where  a  ,a.  ja^, .. .,b  ,b, ,b„, ., .,   are  constants  to  be  determined  so  that  Up 

satisfies   (3.79) •     This  can  be  done,  but  we  shall  not  go  into  the  details. 
We  leave  them  as   an  exercise  for  the  reader. 

The  case  vihere  v  equals  an  integer  is  not  as  easy  to  see  through.     In 
the  next  section  we   shall  discuss  the  singularities  of  a  system  of  n  first- 
order  differential  equations,  and  from  the  results  obtained  there  we   shall 
see  that  vhen  v  equals  an  integer,  the   second  solution  of  (3.77)  must  also 
contain  logarithmic  terms.     VJe   shall  give  further  details    later. 

As  an  illustration  of  an  irregular   singular  point,   consider  the  equa- 
tion 

"         2 
u     +  k  u  =  0 

T  loc  ^i  lor 

at  X  =    00.     The   solutions  of  this  equation  are  u  =  ae         +  pe         ,  which  have 

essential  singularities  at  x  =    oo.     The  singular  point  can  be  transformed 

into  the  origin  by  making  the  change  of  independent  variable  x  =  1/t.     The 

differential  equation  now  becomes 

(3.80)  t^  ^  (t^  ^  )  +  k^u  =  0, 

and  this  has  an  irregular  singular  point  at  t  =  0 .  Note  that  if  we  assume 

u  =  t'^(a  +  a-t  +  ,.,) 
we  will  find  that  there  cannot  exist  values  of  a  ,a  ,a,,  ...  such  that  the 
series  will  satisfy  (3.80). 

To  decide  whether  a  singular  point  x  =  x  of  a  differential  equation  is 
regular  or  irregular  we  proceed  as  follows: 

Substitute  (x-x  )  for  u  in  the  differential  equation  and  in  the  re- 
sulting equation  set  the  coefficient  of  the  lowest  power  of  x  -  x  equal 
to  zero.  The  equation  in  a  so  obtained  is  called  the  indicia!  equation. 
We  state  the  following  Rule,  If  the  degree  of  the  indicial  equation  is 
the  same  as  the  order  of  the  differential  equation^  then  x  =  x^  is  a  re- 
gular singular  point.  If  the  degree  of  the  indicia!  equation  is  less  than 
the  order  of  the  differential  equation,  then  x  °  x^  is  an  irregulrr  singular 
point , 
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We  shall  illustrate  this  rule  by  a  few  examples.     Substitute  u  =  x° 

into  (3.77).      The  lowest  power  of  x  in  the  equation  so  obtained  is  x*     , 

2 
whose  coefficient  is  a(a-l)   +  a  -  v     =0.       Since  the  degree  of  the  in- 

dicial  equation  is  two,  which  is  also  the  order  of  the  differential  equa- 
tion, we  conclude  that  the  point  x  =  0  is  a  regular  singular  point  for 
(3.77).     For  our  second  example  we  take  equation   (3.80).     If  we  put 
u  =•  t     into  (3.80),  we  find  that  the  indicial  equation  is  k     -  0.     If 
k  is  not  zero,  this  equation  is  impossible,   and  therefore  the  point  t  =  0 
is  an  irregular  singular  point  for  (3.8C).     As  our  final  example  we  con- 
sider the  equation 

(x^u')'+  u'   -  0    . 
The  indicial  equation  is  now  a  =  0.     Since   this  is  an  equation  of  first 
degree  while  the  differential  equation  is  of  the   second  order,   we  con- 
clude  again  that  x  =  0  is  an  irregular  singular  point. 

Note  that  to  determine  whether  infinity  is  a  regular  or  irreg-ular 
singular  point,  we  replace  x  by  1/t  and  use  the  rule  to   test  whether  the 
point  zero  is  a  regular  or  irregular  singular  point. 

We  shall  not  give  a  complete  proof  of  this  rule,  but  we  shall  state 

a  theorem  which  partly  proves  the  sufficiency  of  the  rule.     The  proof  of 

this  theorem  will  follow  from  Theorem  III-l  of  Appendix  III. 

Theorem  X.     If  a.   is  a  root  cf  the  indicial  equation  such  that  a  +  ra  is 
o o 

not  a  root  (m  is  any  positive  integer),  then  the  differential  equation 

has  a  solution  of  the  folloxd.ng  form; 

"o   ,  2  , 

u  =  X       (a  +  a^x  +  a^x     +    ...)    . 

Problems ; 

3.30. Expand  the   solutions  of  the  hyper geometric  equation 

x(l-x)  u"  +    f  Y  -  (a+P+1)  x3  u'   -  apu  »  0 
in  povrer  series  about  the  origin. 
3. 31. Expand  the   solutions  of  Legendre's  equation 

[  (l-x^)u']  '  +  vu  -  0 
in  power  series  around  the  point  x  =  1.     Hint.     Put  x  =  1+t  in  the  differ- 

n  0 

ential  equation.     Assume  u  =  t  (1+a^  t+a„t    +    ...),   and  continue  as  for 
Bessel's  equation. 
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3. 32. Consider  the  equation 

u»  +  q-.u'  +  QpU  =  0 
where  q, (x),  q^Cx)  are  each  the  quotients  of  two  power  series  in  x.  Show 

that  the  indicial  equation  at  x  =  0  has  degree  two  if  and  only  if  xq, (x) 

2 
and  X  qp(x)  are  both  bounded  as  x  approaches  zero. 

3.33. "Consider  the  equation  u"  +  q-,u'  +  q^u  =  0.  Suppose  that  in  the 

neighborhood  of  infinity,  q,  and  q^  can  be  expanded  in  power  series  in 

1/x.  WJiat  conditions  on  the  power  series  will  imply  that  infinity  is  a 
regular  singular  point  of  the  differential  equation? 

Illustrative  Example  -  Logarithmic  Singularity 

In  case  the  conditions  of  Theorem  X  are  not  satisfied,  that  is,  if 
X  =  m  is  also  a  root  of  the  indicial  equation,  then  the  solutions  of  the 

differential  equation  may  have  logarithmic  singularities.  The  general 
theory  of  such  cases  is  discussed  in  Appendix  III.  We  shall  illustrate 
this  theory  by  finding  expansions  near  x  =  0  for  the  solutions  of  the 
equation 

(3.81)  u"  +  i  u'  +  (1  -  ^  )  u  =  0  . 

X 

Note  that  this  is  equation  (3.77)  for  the  special  case  v  =  1. 

The  indicial  equation  is 

2 
a(a-l)  +a-l=0  or  a  =1. 

Vfe  have  already  obtained  the   solution  associated  with  a  =  1: 

2  i; 

^1  "  ^^-^  '  2.2-2     *     li.2.^.^.3.2  -...). 

Me  shall  now  find  the  solution  associated  with  a  =  -1. 

In  accordance  with  the  theory  discussed  in  Appendix  III,  we  assume 

—1  2 

u  =  x~  (a  +a.-x+a_x  +...)  +  u,  (x)  log  x 
0  1   4:  i 

■T"  n-1  .  ,     V  1.  n+1 

=  •■  P  X    +  lof  X  /_  >)  X 

* —  n           * —  n 

0  0 


-  155  - 
We  hpve 


00  o  00  , 

n-3  ^    -,  X~  /    .  -I  \  -u     n-1 


Ti"=  7"  (n-l)(n-2)a  x^~-  +   lo^  x  T"  (n+l)nl)  x 
' —  n.  ■*"  n 

o 

+  21  (2n+l)b  x' 


o  o 

n-1 

r 
n 


0 
00  _ 

u«  =  7"  (n-l)a  T^''^  +  log  I  7"  (n+l)l)  x'^  +  7"  b  x*^. 
*—  n  ^~  n  *—     n 

0 

Su'bstitutin^  these  expressions   in  the  differential  eq,us.tion,  we  get 

00  o  -^  -1         00 


T.  (n-l)(n-2)a  x'^'^+logx  ti  (n+l)nb  x"-^+  T  (2n+l)b  x"""^ 

0  0  0 

+  H  (n-l)a  x'^'^+logx  r  (n+l)l)  x°-^+  ^  b  x''"^ 
*—  n  "^     *—  ^        '   n  * —    n 

O  0  0 

v^    n  1     ^    n+1  ^    n  3     r^ 
+  X  »„3f  '  +logx  21  T^„3c   -XI  a„i  "  -logx  2_  ^ 


^    n-r ,    ^  ,   1^1  5?    n-3  ,    ^  ^  n-1 

0  0  0  0 


00 


=  II  ^'^~  &(n-2)a  +a  p  +  (2n-2)b  /] 


n  nr2         n-2- 


05     . 

+  log  X  21  x'^"  [n(n+2)b  +-b„  oH  =  0- 
*"~     ^      n  n— 2 

0 


Since  this  must  "be  identically  zero  for  all  values  of  x,  the  coefficients 
of  X  ~  logx  and  of  x  "  must  be  zero.  This  gives  the  following  set  of 

equntions  for  a  and  b  : 

n      n 

n(n+2)b  +  h    =  0, 
n    n-<d 

n(n-2)a  +a  ^  +  (2n-2)b  „  =  0. 
n   n-2         n-2 


These  equations  are  to  hold  for  n  =  0,1,2,...*  Note  that  we  have  assumed 

^-1  =  ^-2  =  ^-1  =  ""-Z   =  °- 

It  is  clear  that  all  the  odd-numbered  coefficients,  both  aon+i  ^^^ 
b       , ,  are  zero.     We  now  determine  the  even-numbered  coefficients. 
Put  n  =  2  and  we  get 

2  •  k-b^    +     b^     =     0, 

ao     +   2bo  =  0. 
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If  we  essiime  1)  =  1,  then  a  =  -2,  "b-  =  -  :rrr  .  Put  n  =  ^  and  the 
0  o     '   2     2*4 

equat  i  oiB  1)60  ome 

4«6  a^^  +  ag  +  6h2  =  0. 

The  first  equation  determines  t.  ,  "but  the  second  eqviation  determines  neither 

B.^   nor  s.^.      In  fact  we  may  prescribe  a.  arhitrarily,  determine  a^  as 

-(a2+6'b_)/4»2,  and  then  calc^ilate  successively  all  the  other  coefficients  in 

the  expansion  of  u(i). 

This  shows  that  the  expansion  of  u(x)  is  not  unique  since  we  can  always 

add  to  u(x)  an  arhitrary  multiple  of  u^ (x)  ans  still  have  a  solution  of  (3»6l) 

whose  expansion  starts  with  a  term  in  x  .   In  practice,  this  arbitrary  miiltiple 

of  u_ (x)  is  so  determined  that  the  solution  u(x)  will  have  certain  properties, 

for  exajuple  a  given  "behavior  as  x  approaches  infinity.  Meanwhile,  we  assume 

ap  =  0  and  obtain  finally  the  expansion 

k  2 

Up  =  x~-^(-2  +  -~--r -...)  +  z  log  x(l  -  2L_+  ...). 

Problems, 

3.3^.  Find  the  expansion  aro\ind  x  =  0  of  two  linearly  independent  solutions 
of  the  system        .  , 

\ 

U2  /      I  -X 

3.35.  'ind  the  expansion  around  x  =  0  of  that  solution  of 

T  2 

u"+  -u«+(l-S_)ii  =  0 

*  x^ 

which  begins  with  x"'^. 

Solution  at  An  Irregular  Singular  Point 

In  most  physical  applications  of  differential  equations,  and  especially 
in  prcblems  dealin.?  with  wave  propagation,  the  point  at  infinity  is  an  irreg- 
ular singular  point.  As  an  illustration,  consider  the  equation 

u"  +  k^  u  =  0, 


where  k  is  an  arbitrnry  real  number.  The  solutions  of  this  equation  are  of 

+ikx 
the  form  u  =  cos  kx,  u  =  sin  kx,  or  u  «»  e— 

tions  have  an  essential  singulsTity  at  x  =  co  . 


+ikx 
the  form  u  =  cos  kx,  u  =  sin  kx,  or  u  «»  e—    ,  Notice  that  all  these  solu- 
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We  shall  consider  a  more  general  equation, 

(3.82)  u"  +  (k^  +  q(x))  u  =  0, 

and  we  shall  show  that   if  q(x)   is  small  enough  at  infinity,  or  more  pre- 
cisely, if 

OD 

(3.83)  /        |q(x)|djc<00, 

then  thei^  exist  solutions  of  (3.82)  which  behave-  like  e^"^^  at  x  =  oo  . 
This  last  phrase  is  to  be  understood  in  the  following  sense: 

There  exists  a  function  u^(x)  which  is  a  solution  of  (3.82)  and  which 
is  such  that 

lim       |e~       u-(x)  -  l|    =0, 
x->oo 

Similarly,   there  exists  a  function  u„(x)  which  is  a  solution  of  (3.82) 
and  which  is  such  that 

lim       |e      U2(x)  -  l|   =0, 
x->oo 

To  prove  these  results  substitute  u  «  ve"^  in  (3.82)  and  we  get 

v  +  2ikv  +  qv  a  0 . 

If  we  can  now  show  that  this  equation  has  a  solution  v^(x)  which  approaches 

unity  as  x  goes  to  infinity,  we  may  put  vlAx)  =  v^(x)e         and  will  have 

established  the  existence  of  the  function  u^(x).     Similsrly,  by  putting 

••ikx 
u  =  ve  we  can  prove  the  existence  of  u_(x)» 

The  method  of  proof  is  smilar  to  that  used  in  the  proof  of  Theorem  II-I. 

We  write  the  differential  equstion  in  the  form 

It  t 

V     +  2ikv     =  -qv 

and  consider  the  right-hand  side  as  a  non-homogeneous  termj   then  by  using 
the  Green's  function  for  the  left-hand  side  we  shall  obtain  an  integral  equa- 
tion for  v(x). 

The  appropriate  Green's  function  can  be  found  from  Theorem  vni  if  we 
recall  that  v(oo  )  =  1  and  v  (oo)  ■  0.     We  have 


(3.8U)  v(x)  -  1  -  /   i=2__ q(£)v(?)d?. 

CD 

So  far,  these  results  were  obtained  on  the  assumption  that  there  exists  a 
function  u^{x)   and  a  corresponding  function  v(x).  We  shell  shortly  construct 
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a  function  v(x)  which  satisfies  (3.81;).  Then  it  can  be  shown  by  differentiation 
that  v(x)  satisfies  the  differential  equation,  Sinos  it  follows  from  (3.8U)  thst 
y(x)  converges  to  unity  as  x  approaches  infinity,  we  will  have  established  the 
existence  of  u, (x). 

The  existence  of  a  solution  v(x)  of  (3.8U)  will  be  established  by  the  iteration 
method  used  in  proving  Theorem  II -1,  We  put 


,x  ,_-2ik(x-4) 

00 


%+l^''^  °  ^  " «/"   ^"^  2iJc  "  ^'     q(^)V^^'^^   (n  =  1,2,3...) 
Vjix)   =  1, 


and  we  shall  show  that  v  (x)  converges  to  a  limit  function.  We  have 

X  -,  -2ik(x-C) 
V2(x)  =  1  -  /  -"^  2^ q(£)df. 

00 

Note   that  because  of  (3.83)   the  inte^al  convers^es  for  all  real  values  of  kj 

therefore  Vp(x)   is  bounded  for  all  values  of  x,     3y  induction  it  is  easy  to 

prove  in  the  same  way  that  v  (x)   is  bounded  for  all  values  of  x. 

Consider  the  difference  v     t(x)  -  v  (x).     We  have 

n+±      n 

f  ,  -2ik(x-4) 

^^•85)     v^^(x)  =  -  /  i=2___ q(r)[v^(C)-v^_^(C)]dt   n  -  2,3,.., 

CO 


^          -2ik(x-£;) 
V2(x)-v^(x)  -  -  /  -^"^  2ik '^^^^'^^     • 


00 


Since  |exp^2ik(x-^)j  |   =  1  for  all  real  values  of  k,  we  conclude  that 

X       X 

(3.86)       |v  (x)  -V  (x)|  <  /  I  /   |q(C)|  dl 


oo     oo 


Put  X 

Q(x)  =  /   |q(C)|d?  , 

00 

then  (3.86)  becomes 

|v2(x)  -v^(x)|  <^. 

We  now  prove  by  induction  that 

(3.87)  IViW-V==)l::^ 

nlk 
Suppose  this  inequality  holds  for  n  =  l,2,,,,,m-l|  then  from  (3«85)  we  Ret 
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vin-1 


00       (ra-l)lk        (m-Dl  k   oo 


_  Q(xr 

mlk™ 

which  proves  (3.37)» 

Ihe  infinite  series 

(3.88)  v^(x)  +  (y^ix)-v^{x))   +  (v^(x)-V2(x))  +  ... 

converges  for  all  values  of  x  because  each  term  is  less  in  absolute  value 
than  the  corresponding  term  of  the  absolutely  convergent  series 

,  ^  Q(x)  ^  Q(x)^  ^      Q(x)"  ^ 

^     1^  nlk^ 

Since  the  partial  sxms   of  the  series  (3.88)  are  v  (x),  this  proves  that  v  (x) 
■^  n  n 

converges  to  a  limit,  which  we  denote  by  v(x).  It  is  easy  now  to  show  that 
v(x)  satisfies  the  inte^^ral  equation  (3.8U)  and  that  v(oo)  =  1,  As  mentioned 
above,  the  existence  of  the  function  u.  (x)  is  an  immediate  consequence  of  the 
existence  of  v(x) . 

We  state  this  result  as 
Theorem  XI.  If  k  is  real  and  if 

00 

y      |q(x)  I  dx  <   00, 
then  there  exists  a  solution  u,(x)   of  the   differential  equation 

u"*   (k+q)u  =  0 


such  that 

lim      ly-ie"         -I  J 


0. 


X->00 


There  also  exists  a  solution  Up(x)   such  that 

lim     j  Upe'''       -ll      =0. 
x->oo   ^  -^ 

By  making  suitable  transformations  we  can  apply  Theorem  XI  to  much 

wider  dass  of  differential  equations  than  its  appearance  suggests.  For 

exarple,  consider  the  differential  equation 

2    v^ 
(3.89)  (xu')'  +  (k  X  -  -^  )  u  =  0. 
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If  we  put  u  »  x"  '  w  into  (3.89)  we  get 

w"  +  (k^  -  -^LiVk  )  w  =  0. 

X 

From  Theorem  XI  it  is  clear  that  this  equation  has  solutions  which  behave 

-1  lor      ^ikx 
at  infinity  like  linear  combinations  of  e    and  e     ;  consequently  the 

-1/2 
solutions  of  (3.89)  behave  at  infinity  like  x  '   times  linear  combinations 

_  ikx   ,  -ikx 
of  e    and  e 

Problems . 

3.36.     Show  that  vhen  the   substitution  u  =  vexp  L-l/2ybdx  J  is  made  in  the 

differential  equation  u"  +  b(x)  u'   +  c(x)u  =0,  the   resulting  equation  for  v 
does  not   contain  any  first  derivative.     Apply  this   substitution  to  the  equa- 

tion  u''  +  -u'  +ku  =  0,   and  then  by  the  use  of  Theorem  XI  obtain  the  asymp- 
totic behavior  of  u  for  large  values  of  x, 

3.37<     Consider  the  equation  -(pu')     +  aqu  -  pwu  =  0,  where  a  and  p  are  con- 
stants. Fut  v=gu,   g=(wp)'      jy=/     (w/p)         dx,  and  then  show  that 

.iz^   (Q-P)  v  =  0,  where   Q=i    1^     *^   . 
dy^  ^     dy''         ^ 

3.38.  If  /      \Q.\   dy  <oo    ,  where   Q  is    the  function  defined  in  Problem  3.37,  then 
Theorem  XI  may  be  applied  to  the  equation  in  v.     Show  that   in  this  case  we  have 

the   asymptotic  formula  u  =  g~     exp  [+  ip  '     y     (w/p)   '      <i'^  J  • 

2 

3.39.  Find  the  behavior  of  u  for  large  Tallies  of  x  if  -u"  -  px  u  +  au  =  0. 

Hint.     Use  Problem  3.38, 

3.U0       Find  the  behavior  of  u  for  large  values  of  x  if  -u"  =   (x+\)u.     Hint. 

Use  Problem  3.39  with  q  =  0,  w  =  x+X.  . 

Boundary  Conditions  at  a  Sin-^-ilar  FointT 

In  the  previous  sections  the  behavior  of   the   solutions  of  the  homo- 
geneous differential  equation  at  a  singular  point  was  discussed.     V/hat  can 
be  said  about  the   solutions  of  the   corresponding  non-homogeneous  equations? 
Or,   an  equivalent  question,  what  kind  of  boundary  conditions  may  we  reason- 
ably impose  upon  the  solutions  of  non-homogeneous  equations?     That  this  ques- 
tion is  not  trivial  may  be  seen  from  the  folloi-dng  simple  illustrative  problem: 

t  A  similar  discussion  is   given  by  F,  Rellich,  Matheraatishe  Zeitschrift, 
U9,  p.  702  (19UU). 
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Find  a  function  u(x)  such  that  it  is  a  solution  of  the  equation 
u"  +  ^  u'  =  1 

X 

and  such   thst  u(0)   =   1  and  u(l)  =  0. 

We  shall  prove   that   it   is   impossible  to  find  such  a  ftonction  u(x).      The 
solutions  of  the  homo^^eneous   differential  equation  are   1  and  t     .      Since  a 
particular  solution  of  the  non-homogeneous   differential  eq.-uation   is  i  /^, 
the  general  solution  is 

u  =  x^/iJ  +  a  +  hx""^, 

where  a  and  h   are  arbitrary  constants.  To  satisfy  the  boundary  condition 
u(0)  =  1  we  must  have  e  =  1  and  b  =  0,  but  then  u(l)  =  5/^  and  not  zero; 
consequently  it  is  Impossible  to  satisfy  both  boundrry  conditions  simi;l- 
tPneously. 

i^f  course,  the  trouble  is  due  to  the  singularity  at  x  =  0.   However, 
the  followiniF  v^ry  similar  problem  has  a  solution:  find  u(x)  such  that 

u"-  -  u'  =  1 

X 

pnd  snc^  thpt  u(0)  =  1  and  u(l)  =  0.     The  solution  is 

12  2 

■a  =  —  X  '  log,  X  +   1  -  X   . 

We  wish  to  investigate  this  phenomenom  in  order  to  be  able  to  distir^fuish 
between  the  cases  when  a  boundary  condition  can  be  satisfied  and  when  it  can 
not.  For  simplicity  we  shell  treat  only  second-order  differential  operators, 
but  most  of  the  results  can  be  generalized  to  n-th-order  differential  oper- 
ators and  to  systems  of  first-order  differential  equations. 

Any  second-order  differential  equation  may  be  written  in  the  form 

(3.90)        i  [I-(pu')'+q^  =  f(x). 

We  assume  that  p(x),  q(x),  and  w(x)  are  piecewise  analytic  functions  of  x 
in  the  interval  (0,l),  and  we  also  assume  that  p(x)  and  w(x)  are  non-nega- 
tive in  that  interval.  We  suppose  that  the  differential  equation  has  a 
sine;nlarity  et  x  =  0  and  that  hAy)   and  u  (x)  are  solutions  of  the  homo- 
geneous differential  equstio-> 

-  (pu')  '  +  QU  =  0. 
We  suppose  also  thpt  u^(x)  and  uAx)   are  so  normalized  that  the  conjunct 
of  u.  and  u_  is  unity. 

Any  solution  of  the  homogeneous  equation  will  be  expressible  in  the  form 
u  =  au^  +  bu. 
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whTC  a  and  "b  are  ar'hltrsry  constpnts.  We  can  determine  the  values  of  a  and 
"b  "by  considering  the  conjunct  of  m  with  the  fTinctlons  u.  and  u.  in  turn.  We 
hpve 

J(n,u- )  =  sJ(u,Uj^)  +  hJCu  ,U]^)  =  -h 
nnd 

JCn.u^)  =  aJCu^.u^)  +  h  J(n^,ix^)   =  e. 

These  eauptlona  show  thpt  we  msy  define  the  solution  \i  of  the  differen- 
tial eqnstion  hy  Tnepns  of  the  values  of  th**  conjunct  of  u  with  particular 
solutions  of  that  equation.  In  many  casps  it  is  not  necessery 

to  Vnow  the  -oarticular  solutions  completely;  it  is  enough  to  know  the  first 
terms  in  their  exTansions.  For  example,  consider  the  differential  equation 

(3.91)         (xd')'  +  k^  xu  =  0. 

The  point  x  =  0  is  a  regular  singular  point.  By  the  theory  of  the  preceding 
sections  we  know  that  there  is  a  solution  of  the  equp.tion  (3.91)  whose  ex- 
pansion in  the  neighhorhood  of  x  =  0  begins  with  unity,  and  that  there  is 
another  solution  whose  expansion  in  the  neighborhood  of  x  =  0  begins  with 
log  X.  Consequently,  an  appropriate  boundary  condition  at  x  =  0  for  a 
solution  of  (?.9l)  would  be  that 

lim  J(u,a+blogx)  =  0, 
x-^OD 

or,  written  out  in  full,    that 

lim  X   [tt- -   u'(s+blogx3    =  lim (bu  -  fxii '-bxu'lorx|=  0. 
Ei^re  again  a  and  h  are  arbitrary  constants.      If  we  put  b  =   0,    then 


lim     xu'   =  0 
x->  0 

would   >>e   the  appropriate  hotindpry  condition  for  a    solution  of    (3.91 )   that 
is  regular  at  x  =  0. 

It   is   Interesting  to  note   that   the  typical  boundary  condition  at  a 
regular  point,   namely 

ou  +   Pu'  =  0, 
can  be  p^it   in  the   form  of  a   con.Vjnct   condition.      If  x     is   a   regular  uoint, 

0 

then  p(x  )  ^  0  and  there  exist  functions  u.  (x)  and  vlAt)   which  are  solutions 
of  (3.90)  for  f  =  0  and  which  are  also  such  that 
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VL,(X    )    =»    1,    VLiix)    =    0 
1       O  10 

vlAt   )  =  0,  u'(x   )  =  1   , 
Z     0  <c     0 

If  we   tnice  11   (t)  =  -Bn.  (y)  +  oni-Cx),    then  the  typical  botmdpry  condition 

o  l  ti 

given  above   can  be  written  ns   the  following:  conjiinct   condition:  , 

J(u,u  )  =  0. 

0 

Note  again  that  we  do  not  need  to  kaiow  u  (x)  completely, hecaiase  this  condi- 
tion is  equivalent  to  the  following: 

lim   j[u,-P4-a(x-i  D  =  0- 
X— >  X 

0 

Similar  results  hold  in  the  neighborhood  of  an  irregular  singular 
point.  Consider  the  equation 

u"  +  (k^+q)  u  =  0 

and  the  solutions  u^(x)  and  UgCx)  defined  in  Theorem  XI.  Any  solution  u  of 
this  equation  will  be  p  linear  combination  of  u,  end  u_;  consequently  it 
may  be  written  as  n  =  cm.  +  Ptu.  However,  from  the  results  of  Theorem  XI 
it  is  clear  that  we  may  define  u  by  the  two  conjunct  conditions 

lim  J[u,e^^^  =  -2ik0, 
x->oo 

lijn  J&.e'^^I]  =  2ika. 
x->oo 

We  may  e:q)re8s  the  results  of  otir  discussion,  in  a  siniple  form  if  we  intro- 
duce the  concept  of  the  initial  of  a  solution  of  a  differential  equation.  We 
define  the  Initisd  by  considering  the  different  cases.  At  an  ordinary  point  the 
initial  of  a  solution  is  the  constant  and  linear  term  in  the  power  series  expan- 
sion of  the  solution  about  the  ordinary  point.  At  a  regular  singular  point, 
asstDue  that  the  indicial  equation  has  roots  a  =  og.  and  a  =  02*  then  the  differ- 
ential equation  will  have  one  solution  m.   whose  eacpansion  begins  with  x  •*•  and 
other  solutions  n_  whose  expansion  begins  with  x  2.  ^le  general  solution  of  the 


*2 


L^   +  &^^ 


differential  equation  will  he  u  =  P.u.  +  ^^2*    *®  °^^  %^   *  ^2^   *^®  initial 

of  the  solution,  ITote  that  if  a,  =  ou,,  the  expansion  of  u^  will  begin  with 
a  ^    ^       ou, 

X  1  log  X.  In  that  case  the  initial  of  u  is  x  J•(^.+  P^  ^^S  3c).  Theorem  II 

shows  that  an  irregular  singularity, 


-  16if  - 

ikx     — ikx 
there  exist  solutions  ^x^   »nd  Uj,  which  "behsve  like  e    and  e     respectively. 

The  general  solution  will  he  of  the  form  u  =  P^u^+P_u  .  V/e  call  P  e    + 

-ikx  i  i  <i  2  1 

^2®     the  initial  of  the  solution.  Note  that  according  to  the  discussion 

of  the  two  previous  sections,  the  initial  of  the  eolution  majr  "be  fonmd  quite 

easily^  even  without  knowledge  of  the  complete  solution.   Now,  we  state  the 

Rule.  A  B-'j.ita'hle  "boundary  condition  at  y  =  y_  for  the  solution  u  of  a 

homogeneo'is  differential  equation  is  ohtained  hy  giving  the  value  at  x  =  x 
— ^ — ^ —      . __^^___  o 

of  the  con.lunct  of  u  with  the  initial  of  any  particular  solution  of  that. 

equation. 

^ound!"ry  Conditions  for  Non-Homogeneo-'JB  Differential  Equations 

In  this  section  we  shall  discuss  the  possihle  boundary  conditions  that 
may  he  inroosed  on  s  non-homogeneons  differential  equation.  The  essential 
reason  for  the  distinction  between  homogeneous  and  non-homogeneous  equations 
is  thst  in  the  homogeneous  case  it  is  immaterial  how  sin^rular  the  solutions 
are,  while  in  the  non-homogeneous  cafie  the  solutions  are  defined  "by   integrals 
which  may  not  converge  at  a  singular  point.  Because  of  this,  we  shall  have 
to  distingiiish  "between  the  cases  where  the  solutions  of  the  differential 
equation  are  intpg-^a"ble  and  the  cas^s  where  the  solutions  are  not  integrahle. 

Consider  again  the  equation  (3.90): 

I«i  =  7  C-  (pu')'  +  quU  =  fix) 

in  the  neighborhood  of  x  =  0.  We  assume  first  that  i  =  0  is  an  ordinary 
point  of  the  equation.  This  Implies  that  p,  q,  and  w  are  regular  at  X  =  0 
and  that  p(o)  /  0.   In  this  case  it  is  well  known  that  we  may  assign  ar- 
bitrarily the  values  of  u(0)  and  u*(o),  or,  what  is  equivalent,  assign 
the  vrlue  at  x  =  0  of  the  conjunct  of  u  with  the  initials  of  any  two  solu- 
tions of  the  homogeneous  equation. 

Supijose  now  that  x  =  0  Is  a  singular  point,  and  let  m-  (x)  and  u„(x) 
be  independent  solutions  of 

-  (pu')  •  +  qu  =  0. 

We  assume  that  JCu-.Ug)  =  !•  '^^^     Green's  function  for  the  operator  L  on 
the  manifold  where  u(o)  =  u'(o)  =  0  is 

Ei^CxH^C!)  -  u^U^Clfj]  H(x-  f ); 
hence  a  particular  solution  of  (3.90)  is  given  by  the  formula 
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(''.92)        u^  =  .i^(x)/u2(t)w(f)f(f)df   -U2(x)/\(()w(?)f(f)df  . 

By  dl-PferentiatiDn  we  find 

u^  =  tL[(T)/u2(^)w(hf(f)df   -  u^(x)/\(f^(nf(|)df  . 

yrom  these  formulas  It  Is  easy  to  see  that 

(-^.9^)        J(u^,u^)  =  -J(m^,m^)   ^^ix)  =   P^^x). 

J(u  ,u  )  =  J(n^,u  )  Pj^(x)  =  Pj^d), 

where 

P.U)  =/u.(]p)w(!F)f(!F)df  . 
-i     1  -^ 

P^(x)  =/u2(f)w()F)f(f)  df  . 

We  shall  now  show  that  if  these  integrals  converge  as  x  approaches 
zero,  we  may  impose  any  two  conjunct  conditions  on  u(x)  at  x  =  0.  For 
example,  suppose  that  the  conditions  on  the  solution  of  (?.90)  are 

J(u,  a^,a^  +  Va^^o  ^  °1  • 

J(u,  a  u,  +  h-u.)    =  c„. 
^  ^    '^  ^  x=0    ^ 

Any  solution  of   ('^.90)  may  he  written   in  the   form 

^  =  ^  ■"  ""l^  *  ''2^2   • 

where  V     and  *'     are  arhitrary  constants.     With  the  help  of   (3.9'')  we  find 
that 

J(V^V^2^2'^lV^^2^^=  ^^2  -^  \h  *  "^1^2  -  '^2^  =  '^l' 
'^K*^VV2'«2V^2^^^0  =  «^2P2"*-^2Pl  *  V2  -  V2  '  =2   » 

where  we  have  put  P^  =  3^(0),  P^  =  ^2^°^* 

When  these  line,-^r  equations  are  solved  for  V  and  V   the  formula 
u  =  '^■q'*''^ i'\'*''^ 2^2  ^^"^^^  ^   function  u(x)  which  satisfies  the  differential 
eq\3ation  (3»90)  and  also  the  given  conjunct  conditions.  The  equations  for 

"^1  and  V2  will  have  a  unique  solution  if  and  only  if  a^bp  _  a^h  ^  0. 
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This  condition  merely  expresses  the  requirement  that  the  functions 
a,u^+b-U2  and  a^iLj+bpUp  be  linearly  independent. 

Note  again  that  we   may  VTite  the  conjunct  conditions  in  terms  of  the 
initial  of  the  solutions  instead  of  the  full  solutions. 

We  must  now  investigate  in  what  cases  the  integrals  in  (3»92)  will 

converge.  One  special  case  is  of  particular  importance.  Consider  the 

scalar  product  with  w(x)  as  a  weight  factor,  that  is, 

1 
(u,v)  =  /*  u(x)v(x)w(x)dx. 

Note  that  with  this  scalar  product,  Lu  is  formally  self-adjoint  and  the 
integrals  up  to  x  =  0  in  (3.92)   are   (Up,f)   and  (u,,f).     We   say  that  a 
function  u  belongs  tOoCp  if  (u,uj  is  finite.     If  u  and  v  both  belong 
to  (Ap  then  (u,v)   is  finite  because  of  Schwartz's  inequality. 
Using  this  resvili  we  conclude  that  if  u,  ,Up,   and  f  belong  to  «>C,,  then  the 
integrals  in  (3.92)   exist  and  any  two  conjunct  conditions  can  oe  imposed 
on  u. 

If  u,    and  Up  belong  to  ofp*    then  since   any  other   solution  of  the 
homogeneous  equation  is  a  linear  combination  of  u^    and  Up,   that   solution 
also  belongs  too<p.     We   shall  say  that  the  operator  L  is  of  the  limit- 
circle'     type  to  X  =  0  if  all  the  solutions  of  Lu  =  0  belong  to  ^p  at 
x  =  0.     We  now  state  the  result  of  our  discussion  in  the  following 
Theorem  XII .     Sup-pose   that  the  operator  L  is  of  the  limit-circle  type  at 
X  =  0,   and  suppose  that  f(x)  belongs  to  y^,  at  x  =  0;   then  there  exists 

a  function  u(x)   in  <JC'n  which  satisfies  the  equation 

»i 

Lu  =   f 

and  which  satisfies  any  two  independent  conjunct  conditions  at  x  »  0. 

As  an  illustration,   consider  the  operator 

Lu  =  — :  (xu' )'    -  k  u. 

If  1 

\{f{x))^  dx  <  oo. 


/■■ 


there  exists  a  solution  u(x)  of  the  differential  equation 
-Lu  =  u«'  +  -  u'   +  k^u  =  f 

X 

which   satisfies  the  two  conditions: 

1"  The  name  and  the  concept  were  first  introduced  by  H.  Weyl,    "  tJber  geWohnliche 

differential  gleichungen  mit  Singularitaten  und  die  zup;ehorigen  Entwicklungen 
willkiirlicher  Funktionen,"  Math,  y\nnalen,  ^,   220-269  (1910). 
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I 

lim  xu   =  c. 

X  -♦  0 

lim  [u  -  11  X  log  x]  =  Cp. 
x->  0 

Suppose,  however,  that  L  is  not  of  the  limit-circle  type;  then  not  every  solution 

of  La  =  0  belongs  to  ilf  o.  We  shall  then  say  that  L  is  of  the  limit-point''   type 

II   T  I     1  /  3  'v ' 
at  X  =  0.  For  example  if  Ini  =  -u  -  "  u  =  -  — r  ^x  u  ;  ,  then  the  solutions  of 

—2  X 

Irti  =  0  are  u.  =  1  and  u-  =  x  .  Since 

(u,u)  =  y  n  x'dx, 

we  see  that  u^  does  not  belong  to  «C-,  and  that  L  is  of  the  limit-point  type. 

Another  very  important  example  of  a  limit-point  operator,  this  time  at 
X  =  00,  is  the  following: 
(3.9^)  Im  =  -u  -  (k  +  q)ti. 

Fote  that  here  no  solution  of  Iai  =  0  belongs  to  jCo  over  the  infinite  interval. 
We  shall  state  a  theorem  which  will  enable  us  to  discuss  the  appropriate  boundary 
conditions  for  an  operator  such  as  (3.9^) i 

J30 

Theorem  XIII.      If  k  is  real  and  /*°°lq|    dx  <  oo  andy      |f  |dx  <  oo,   then  the  equation 

(3.95)  u"+    (k^  +   q)u     =     f 

has  one  and  only  one  solution  which  approaches  zero  as  x  approaches  infinity. 

The  proof  is  simple.  let  •u.  and  Uo  ^^  t^©  functions  defined  in  Theorem  XV, 

iifimely,  solutions  of 

"   /  2    \ 

u  +  (k  +  q)u  =  0 

liOC     — ikx 
which  behave  like  e    and  e     respectively  for  large  values  of  x.  Note  that 

J(\i^,u  )  =  -2ik  . 

If  we  use  the  Green's  functions  for  the  boundary  conditions  u(oo)=u(oo)=0, 
we  find  that  oo 

^=  /  [u^(x)u2(J)-U2U)u^(f)Jf(i)df  /2ik. 

ikx     —ikx 
Since  u.  and  Up  behave  like  e    and  e     for  large  values  of  x,  they  are  bound- 
ed for  all  values  of  x;  consequently,  u^(x)up(5  )-U2(x)u^(5)  is  bounded  for  all 
values  of  x  and  f  .  We  have  then 

-t^  Also  this  name  ajid  concept  were  first  introduced  "by  H.  Weyl,  ibid. 
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z 
where  C  is  a  3uita"ble  constant.  Since  f  is  integrable,  then  as  x  approaches 
infinity,  the  integral 

^   |f(t)|d 

approaches  zero  and  therefore  u(x)  also  approaches  zero. 

The  function  Viix)   is  a  unique  solution  of  the  pro'blem  hecause  any  other 
solution  of  (3.95)  would  be  of  the  form  u  +  a.u-+a2U2»  where  ol  and  a.  are 
arbitrary  constants,  and  it  is  ohvioxis  that  such  a  solution  does  not  approach 
zero  as  x  goes  to  infinity. 

Note  that  the  theorem  does  not  require  the  function  u(x)  to  he  in  ^  , 
In  fact  this  requirement  cannot  "be  satisfied,  as  may  "be  seen  from  the  equation 

..   2     e^^ 


u  +  k  u  = 


Jlz- 


It  is  easy  to  show  that  for  large  values  of  x,  u  behaves  like  the  function 
ix~  ^  e   /k,  which  is  not  in  oC  2« 

As  an  illustration  of  Theorem  XIII,  consider  the  equation 

(3.96)  u"+  k^u  =  e"^. 

The  only  solution  of  this  eq'oation  which  vanishes  for  large  x  is  u  =  e  (k  +l) 
We  see  then  that  it  is  impossible  to  find  a  solution  of  (3*96)  which  vanishes  at 
infinity  and  also  satisfies  a  boundary  condition  such  as  u(o)  =  0. 

Suppose,  however,  that  k  is  not  real,  but  has  a  small  positive  imaginary 
part.   In  that  case  e   ,  which  is  a  solution  of  the  homogeneous  equation  cor- 
responding to  (3»96),  vanishes  exponentially  as  x  becomes  infinite;  therefore, 
for  any  value  of  a  the  function 

(3.97)  u(x)  =  -f^  +  ae^^^ 

k>l 

is  a  solution  of  (3.96)  which  vanishes  at  infinity.  Now,  we  may  impose  a 
boundary  condition  such  as  ■a(o)  =  0;  then  we  would  find  that  a  =  -(k  +l)~  . 

We  may  clarify  this  result  by  noting  that  when  k  is  real  no  solution 
of  the  homogeneous  equation  corresponding  to  (3.96)  belongs  to  Sfg*  -^^^ 
when  k  has  an  imaginary  part,  one  and  only  one  solution  of  the  homogeneous 


ikx     -ikx 
3    or  e    , 

on  whether  the  imaginary  part  of  k  is  positive  or  negative.     Consequently, 


eqiiation  belongs  to      ^^r     this  solution  is  either  e         or  e  ,   depending 
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when  k  la  real,  the  solution  of  (3.96)  which  vanishes  at  infinity  is  unique- 
ly determined,  and  it  is  impossible  to  impose  any  other  conditions.  But  when 
k  has  an  imaginary  part,  an  arbitrary  multiple  of  e    or  e~    can  be  added 
80  as  to  satisfy  another  boundary  condition. 

In  many  practical  applications  where  an   equation  similar  to  (3,96)  occurs, 
it  is  physically  very  reasonable  to  assume  that  k  has  a  small  imaginary  part. 
This  assumption  is  equivalent  to  taking  into  account  the  effect  of  a  friction- 
al  force  in  mechanics,  finite  conductivity  in  electromagnetic  problems,  or 
viscosity  or  damping  terms  in  acoustic  problems.  Therefore  physical  problems 
which  impose  a  boundary  condition  on  the  solution  of  (3.96)  may  be  solved  by 
assuming  k  to  have  a  small  imaginary  part.   It  is  obviously  immaterial  whether 
k  has  a  positive  or  negative  imaginary  part,  but  for  the  sake  of  convenience 
we  shall  always  asssume  k  has  a  positive  imaginary  part. 

It  is  often  desirable  to  obtain  the  physically  significant  solution  (3.9?) 
without  making  the  assumption  that  k  has  an  imaginary  part.  In  such  a  case  we 
can  no  longer  require  that  the  solution  of  (3.96)  vanish  at  infinity.   Instead 
we  shall  say  that  u  must  behave  like  en  outgoing  wave  at  infinity.  Here 
we  mean  by  "outgoing  wave"  tliat  for  k  real,  u  is  the  limit  of  solutions 
of  (3.96)  which  vanish  at  infinity  when  k  has  a  positive  imaginary  part. 

We  summarize  oijir  discussion  in 
Theorem  XIT.   If  f(x)  and  g(x)  are  integrable  from  0  to  00.  then,  for  real  k. 
the  equation 

(3.98)  u"+  (k^+  q)u  =  f 

has  a  solution  which  behaves  like  an  outgoing  wave  at  infinity  and  satisfies 

any  suitable  boundary  condition  at  x  =  0. 

The  proof  is  based  on  Theorem  XIII  which  states  that  there  exists  a  unique 

solution    u  of  (3.98)  which  vanishes  at  infinity.  Put 
0 

U   =   U   +  CM., 

Where  a  is  an  arbitrary  constant  and  u.  is  defined  in  Theorem  XI  as  that  solution  oi 

u   +  (k  ■♦■  q)u  =  0 

ikx 

which  behaves  liJce  e    for  large  x.  The  value  of  a  can  now  be  determined  so 

th©,t  u  will  satisfy  a  suitable  boundary  condition  at  x  =  0,  This  proves  Theorem 
XIV. 
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Problems ; 

3.^1     Investigate  suitable  1)01011(18X7  conditions  at  x  =  0  for  the  solution 

of 

•'      n     •       ,  2  .,   X 

u    +~ii    +ku    =    f  vx;, 

3»hZ     Consider  the  equation  in  Problem  3.^1  and  discuss  the  'behayior  of 
u  for  large  values  of  x.  Hint.  Remove  the  first-derivative  term  as  in 
Protlem  3.36. 
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Appendix  I 

Green's  Function  for  a  System 

We  muBt  realize  that  since  the  Green's  function  is  the  inverse  operator 
end  since  the  operators  for  systems  are  not  only  differential  operators  but 
also  matrices,  we  must  teJce  for  the  Green's  function  not  a  vector  solution  of 
the  system,  "but  instead  a  matrix  solution.  Consider  the  system  of  n  first- order 
differential  equations: 

with  the  scalar  boundary  conditions 

(1.2)  B-(U)  =  B_(U)  =  ...  B  (U)  =  0. 

The  Green's  matrix  for  this  system  is  defined  to  be  the  matrix  G(x,t)  sxich  that 

(I.?)        s-  ^^  =  **(^-  *^^' 

where  I  is  the  identity  matrix  in  n  dimensions.  Also,  each  column  of  G(x,t)  must 
satisfy  the  boundary  conditions  (1.2). 

It  is  easy  to  see  thst  (1.3)  is  the  definition  suitable  for  solving  the  non- 
homogeneous  equation 
(I.i^)  f2  .  AU  =  P. 


where  P  is  a  given  vector.  We  find  that  the  vector  U  defined  by  the  formula 

U  =  /   G(x,t)p(t)dt 
•Jq 

is  a  solution  of  (l.'O  which  satisfies  the  boundary  conditions  (1,2). 

For  all  values  of  x  j^  t,  the  function  G(x,t)  satisfies  the  homogeneous  equa- 
tion (l.l).  If  we  integrate  (I,?)  with  respect  to  x,  we  get 

G  - 

where  H(x)  is  the  Heaviside  \anit  function.  This  equation  implies  that  the 
diagonal  elements  of  G  have  Junips  of  magnitude  1  at  x  =  t,  while  the  non- 
diagonal  elements  of  G  are  continuous  at  x  =  t.  This  result  can  be  verified 

by  putting 

G  =  +H(x-t)l  +  K(x,t), 

where  K(x,t)  is  a  metrix,  and  substituting  this  expression  in  (1,3).  We  get 

(1.5)  g  -  AZ  =  M, 

where  M  is  a  matrix  with  piecewise  contintious  elements.  It  can  be  shown  that 
(l,5)  has  a  solution  in  which  the  elements  of  K  are  all  continuous  functions,  and 
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conseqnently  the  only  discontinuities  of  G  must  he   those  due  to  the  term 
H(x-t)l.   Since  this  latter  matrix  has  jumps  of  magnitude  1  along  the  diagonal 
when  X  =  t,  the  discontinuity  conditions  on  &  are  verified. 

We  shall  now  give  b  formula  for  the  Green's  matrix  of  a  system.  The  formula 
can  "be  easily  transformed  to  give  the  Green's  function  of  an  n-th-order  differential 
equation. 

Theorem  T-I.   Suppose  that  U-(x),'**,U  (x)  are  n  linearly  Independent  vector 
solntlnns  of  the  system 

(1.1)  f?=A^. 

dx 

and  let  S(i)  he  the  matrix  whose  columns  are  the  vectors  U,  ,♦••, U  .   If  there 
IJ L_n 

does  not  exist  a  solution  H  of  (l.l)  which  also  satisfies  the  hoxindary  conditions 
(1.2),  then  there  exists  a  matrix  C,  whose  elements  are  functions  of  t,  such  that 
the  matrix 

(1.6)         G  =  S(x)C(t)  +  S(x)S(tr-^  H(x-t) 

is  the  Green's  matrix  of  the  differential  equation  (l.l)  with  the  boundary 
conditions  (I.2). 

Note  that  from  the  corollary  to  Theorem  IV  it  follows  that  det  S  }^  0  and  so 
the  inverse,  S(t)~  ,  exists. 

First,  it  is  cl^ar  from  the  definition  (I. 6)  that  each  column  of  G  is  a  linear 

comhination  of  the  columns  of  Jf    •  that  is,  of  U^  ,  •  •  ♦  ,TJ  ,  and  therefore  each  column 

^  in 

of  G  satisfies  (l.l)  for  x  /  t.   Second,  the  change  in  G  as  x  goes  across  the 
value  X  =  t  is  S(t)S(t)~  =  I.  This  means  that  the  diagonal  elements  in  G  have 
Jtunps  of  magnitude  1,  while  the  non-diagonal  elements  are  Continuous.  These 
are  the  correct  conditions  for  the  Green's  matrix  et  x  =  t. 

'P'inally,  we  must  show  thpt  C  cpn  he  determined  so  that  the  boundary  conditions 

(1.2)  pre  satisfied  hy  epch  column  of  G.   Let  CT(t),*",  c  (t)  he  the  elements  of 

i         n 

the  first  column  of  C;  then  the  first  colujin  of  G  will  he 

where  V  represents  the  first  column  of  S(x)S(t)~  .  Since  the  hoTindary  conditions 
are  linear  and  homogeneous,  it  follows  that 

\^^^1^  "  °l\^^l^**"*°n\^^n^'^\^V^^"*^^'  5^1.2. ....n. 

These  n  linear  eauations  for  the  un3cnowns   c,,Cof***iC     will  have  a  solution  if 

i     c  n 

the  determinant  of  the  coefficients,  namely 
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det  Bj^(U^)   (k,j  =  l,2,...,n), 

is  not  zero.  However,  if  the  determinant  were  zero  there  would  exist  constcnts 
a  ,  Op,,,,, a  such  that 

•    0  =  0^3,  (UT)+,,.+a  B,  (U  )-B,  (a  UT  +  .,.+a  U  )  k  -  1,2, ...n. 
1K±  nlcnkil  nn 

This  would  contradict  the  hypothesis  of  Theorem  I-I  because  the  vector 

u^OtU-   +  ...  +a  U 
11  n  n 

would  be   a  solution  of  (I.l)   satisfying  (1.2),       Consequently,   the  determinant 

^     is  not  zero  and  the  quantities   c^,Cp,,,,,c     can  be  found.       Similarly,  we  can 

find  the  elements  of  every  other  column  of  C  and  so  prove  the  theorem. 

It  is  easy  to  show,  either  by  the  use  of  Theorem  I-I  or  by  a  direct  proof, 

that  the  Green's  function  for  an  n-th  order  differential  equation  is  given  by 

the  following: 

Theorem  I-II,     Suppose  that  -il.(x),,,.,p  (x)   are  n  linearly  independent  solutions 

of  ttie  homogeneous  differential  equation 


(1.7) 


T  d  u  ^        ,  p, 

Lu  =  p^    — -  +,,.+  p  u  =  0, 
o      ,  n  -^n  ' 

dx 


If  there  does  not  exist  a  solution  of  (1,7)  which  also  satisfies  the  boundary 
conditions 

(I,3)  B^(u)  =  •••  «  B^(u)  =  0, 

then  there  exist  n  functions  a,  (t), .,.,c  (t)  such  that  the  function 

g(x,t)  =  CtUt(x)  +  .,.C  U  (x) .   QfXj"^)         Hfv-+^ 

11''     n  n^"^'   p^(t)w(u^,'««,u  )  ^'^^'^h 


where 


G(x,t) 


VL^W 

Vi^it) 

• 

• 

• 

\ 


(x) 


Vl, 


M 


u„(t) 


U 


U 


(n-2) 


it) 


n 


is  the  Green's  function  for  equation  (I.7)   and  boundary  conditions  (1,8), 

Here  V7(u^,,,,,u^)  is  the  value  of  the  Wronskian  of  u^(x),.,,,u  (x)  at  x  =  t. 
The  proof  of  this  theorem  is  left  to  the  reader  as  an  exercise. 
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Problems ; 

1,1,       Find  the  Green's  matrix  for  Equation  (1,1)  when 

A  - 


Vk  o  / 


k 

and  the  boundary  conditions  are  that  the  first  component  of  U  equals  0  at 
X  »  0  and  that  the  second  component  of  U  equals  0  at  x  »  1, 

1.2,  Find  the  Green's  function  for  the  equation  u         +k  u  »  0  in  the  space 
defined  by  the  boundary  conditions  u(0)  =  u  (O)  =  u(l)  »  u  (l)  =  0, 

1.3,  Prove  Theorem  I-II, 

I.li,  Let  S  be  the  matrix  defined  in  Theorem  I-I,  Show  that 

dx 
and  that 

^  =  -S"^A 
dx     -^   A, 

Hint,  Differentiate  the  identity  S(x)S(x)~  =  I  with  respect  to  x. 

Appendix  II 

Existence  of  a  Solution  of  Equation  (3»29) 

We  now  wish  to  justify  the  statement  that  the  solution  of  equation  (3«29) 
is  a  contij^aously  differentiable  function  of  x.  We  shall  prove  a  theorem  which 
will  contain  this  result  as  a  special  case. 

Theorem  II-I.  Let  p(x),  q(x),  and  f(x)  be  piecewise  continuous  functions  of  x 
in  the  closed  interval  (0,l)j  and  let  p(x)  be  positive  in  that  interval;  then 
there  exists  a  continuous  function  u(x)  such  that  pu  exists  and  is  continuous 
for  all  xj  furthermore 

(n-l)  (pu')'  -  qu  =  f(x) 

for  all  vaL  ues  of  x  for  which  both  sides  are  continuous  functions  of  x} 

and  u(0)  °  u  (O)  °  0, 

The  proof  will  ^^-^strate  how  the     Green's  function  can  be  used  to  transform 
a  differential  equation  into  an  integral  equation,     (II, l)  may  be  written  in  the  form 

Lu  =  (pu   )     »  qu  +  f(x)  =»  h(x). 
Suppose  that  h(x)  were  knownj  then  u  would  be  found  by  using  the  Green's  function 
for  L,  just  as  the  solution  of  (3.36)  was  ^iven  by  (3,38), 

Consider  the  equation  defining  the  Green's  function,   namely 


-  175  - 


(II.2)  (pg*)'  =  -6(x-t), 


with  g(0)  =  g  (o)  =  0,  It  can  "be  solved  "by  direct  integration.  We  have 

pg   =  -^(x-t)  +  c 
where  c,  the  constant  of  integration,  must  "be  zero*  since  g  (o)  =  0.  Put 

r(z)  =  J      ds/pCs); 

0 

a  second  integration  will  then  give 

g(x)  =  [r(t)-r(x)]H(x-t). 

Multiply  (II. 2)  hy  -h(t)  and  integrate  from  0  to  1.  We  find  that  the  solution 
of  La  =  h  is  ^ 


(11. 3)  ^  =  y  Cr(x)-r(t)]h(t)dt. 

o 
Since  h(t)  =  g(t)u(tKf(t)  is  not  a  known  function  of  t,  equation  (3.84)  be- 
comes the  following  integral  equation  for  u(x): 

X  x 

(II.'O  u(x)  =  J       f(t)[r(x)-r(t)]dt  +  /  q(t)u(t)[r(x)-r(t)Jdt. 

0  0 

We  have  thus  shown  that  if  (II. l)  has  a  solution  u(x),  this  solution  will 
satisfy  the  integral  equation  (II.4),  We  now  show  the  converse.  Suppose  that 

(11.4)  has  a  solution  u(x).  Note  that  from  its  definition,  r(x)  is  a  continuous 

function;  consequently  the  right-hand  side  of  (II. 4),  and  therefore  also  u(x), 

is  a  continuous  function  of  i.  More,  however,  is  time.  When  we  differentiate 

the  right-hand  side  of  (ll.il)  we  find  that 

,        ^x  X 

u  (x)  =  y  f(t)dt/p(x)  +  f    q(t)u(t)dt/p(x). 

0  0 

This  equation  shows  that  for  all  x  such  that  p(x)  is  continuous,  u  (x)  exists  and 

I  ^X  -X 

(11.5)  p(x)u  (x)  =  J     f (t)dt  +  J       q(t)u(t)dt. 

O  0 

Since  the  right-hand  side  of  this  eq\aation  is  a  continuous  function  of  x  for 
all  X,  then  at  the  points  where  p(x)  is  discontinuous  we  may  define  u(x)  so 
that  p(x)u(x)  will  be  continuous  everywhere. 

Because  the  integrands  in  (II, 5)  are  piecewise  continuous  functions  of  x 
we  may  differentiate  the  right-hand  side  of  (II. 5);  we  find  that 

(pu')'  =  f  (x)  +  q(x)u(x) 
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whenever  the  right-hand  side  is  contimiou.8.  This  shows  that u  satisfies  (3.39). 
From  (11,5)  it  is  clear  that  u  (O)  =  0.  This  shows  that  (ll.^)  is  completely 
equivalent  to  (ll.l)  and  its  "boundairy  conditions. 

In  order  to  complete  the  proof  of  Theorem  XII  we  need  only  show  that  there 
exists  a  solution  of  equation  (11,^0.  At  first  glance,  this  does  not  seem  to  be 
a  simplification  since  equation  (II, i|)  seems  more  complicated  than  the  original 
equation  (II. l).  However,  equation  (II. ^)  is  acttially  much  better  suited  for 
our  purpose  because  it  contains  an  integral  operator,  which  is  bounded^ 
while  (II. l)  contains  a  differential  operator,  which  is  always  ^anbounded. 

The  technique  which  we  shall  use  to  prove  the  existence  of  a  solution  of 
(II, 4)  is  exemplifipd  in  the  following  case. 

Suppose  K  is  a  bounded  operator  in  a  complete  linear  space,  and  b  is  a  known 
vector;  we  wish  to  find  a  vector  y  such  that 

(II. 6)  y  =  b  +  ky. 

If  the  bound  m  of  K  is  less  than  one,  we  shall  prove  that  such  a  vector  exists, 
■Replace  (II. 6)  by  the  following  sequence  of  approximating  equations: 

^1   =  ^ 


yirfl=  ^*^n'        ^=  ^'2.,.. 


We  have 


Vl-^n=  ^^^n-^l^ 


for  all  n  >  1,  and  therefore 


Vl-^n    ="'    ^n-Vll-"'^    ^r^-l-^n-zf 

Since 

^P "  ^  '^  ^^P  "Vi^  *  Vi '  V2^*  •••  *^v 

we  have 

^p  "  \^  -  '^p  "  Vi'*' Vi  "  V2  •"  —  *  V: 

m'^lbl 


-V 


<  (mP  +  mP*^+  .,.  +  m*l-^)lbl  =  ^sl  [^1. 


If  m  <  1,  then  m^  and  m'^  approach  zero  as  p  and  q  approach  infinity;  consequently, 
the  limit  of  |y  -  y  I  is  zero.  Since  the  space  is  complete,  the  sequence  y  must 
approach  a  limit  which  we  denote  by  y. 

If  we  let  n  approach  infinity  in  the  equation 
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we  find  that 

y  =  b  +  Ky; 

consequently,  we  hsvc  estaTillshed  the  fact  that  there  exists  a  solution  of 
the  originsl  equation,  namely  the  vector  y,  which  is  the  limit  of  the  sequence 

Now  let  us  return  to  the  proof  of  the  existence  of  a  solution  of  equation 
(11.^).  The  function  u(x)  will  correspond  to  the  vector  y,  the  first  integral 
on  the  right  will  correspond  to  the  vector  h,  and  the  integral  operator  whose 
kernel  is  q(t )[r(x)-r(t^  will  correspond  to  the  hounded  operator  K.  We  proceed 
just  as  in  the  case  of  (ll.6)t 

Put 

u^(x)  =/f(t)Cr(x)-r(t}]dt. 

0 

and  for  n  =  1,2,..., put 

^  ^.M   =  /'f(t)G-(x)-r(tl]dt  +/q(t)[Ir(^)-r(tl]u^(t)dt. 
n+l     </  I'  n 

0  o 

We  have 


vl^At)   -u  (x)  =/q(t)rr(x)-r(tani  (t)-u  .(t^dt 
n+l      n     >^  n     n— ± 

end  therefore 

In^,(x)-u  (x)|  ^/lq(t)I  !r(x)-r(t)|  lu  (t)-u  ,(t)|dt. 
nr  ±  n      <'  n     n— x 


z 

f 

o 

X 

f 

0 

In  order  to  estimate  the  right-hand  side  we  consider  the  function 
(11.8)    Tl(x)  =/:r(x)-r(ta  dt  =/dt/^  =  /^ds, 

O  O     t  0 

where  we  have  used  the  previously  given  definition  of  r(x)  namely 

X 

^^^^  =/  F^   • 

o  ^ 

Since  p(s)  is  positive,  it  follows  that  both  r(x)  and  R(x)  are  non-negative 
and  increasing  functions  of  x.   Note  that 

/R(t)'^Er(x)-r(tadt  =/dtR(t)'^/  ^=/4fy/dtR(t)^ 
0  0        too 

Since  P. (t)  is  ai  increasing  function  of  x. 
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s 


/dtRCt)^  <  R(s)''  /dt  -  sR(s)^  ♦, 
consequently,  if  we  use  the  fact  that  from  (II. 8),  R'(x)  ■  x/p(x),  we  have 

(II.9)y'R(t)^[r(x)-r(t)]  dt  <f  |*y  R(s)^  -y  R'(s)dsR(s)^  =    ^^^  . 

We  can  now  get  an  estimate  of   |u  ^-(x)-u  (x) j  from  (ll.?).     First,  we 

note  that  since  p(x) ,  q(x),   and  f(x)  are  piecewise  continuous  functions  of  x, 
they  must  be  bovinded  in  absolute  value;     also,   since  p(x)  is  positive  in  the 
interval  (0,1)  it  is  greater  than  some  fixed  constant,   and  consequently  r(x) 
is  bounded.     Let  m  be  a  bound  for  the  functions  p(x),  q(x),f(x),  and  r(x); 
then  from  the  definition  we  have 

lu  (x)    I  </  |f(t)|    [r(x)-r(t)  ]  dt  <  m  R(x), 

^0 

and  from  (II. 7)  and  (II.9)  we  find  that 

|u2(x)-u^(x)  I  <f  |q(t)|    [r(x)-r(t)  ]mR(t)dt 
o 

<  m2y'R(t)[r(x)-r(t)  ]  dt  <  ^J^  . 
o 

Now  we  use  mathematical  induction.     Suppose  that  for  n  «  l,2,...,k  , 

|u   (x)-u     ,(x)    I  <     !l^^l!L  ; 
'  n  n-1  '  nl  ' 

we  shall  prove  the  same  result  holds  fcr  n  =  k+1.     From  (II. 7)  and  (II. 9)  we  have 
|Uj^^-l(x)-Uj^(x)|  </     |q(t)l    [r(x)-r(t)]  m^R(t)*'  dt/ki 


<E!;^/H(t)^[r(x)-r(t)]dt<2!:^j} 


0 

Consider  the  infinite  series 

u^(x)  +  [  U2(x)-u-j^(x)  ]  +  [u^(x)  U2(x)  ]  +  ... 

We  have  just  shown  that  each  terrn  of  this  series  is  less  in  absolute  value  than 
the  corresponding  term  of  the  series 
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mE(x)  +  ^^Ij^  +  S^^  +  ...  =  -1  +  e3p[niE(x)]. 

This  latter  series  converges  absolutely  and  uniformly  for  all  values  of  x  in  (0,1); 
consequently  the  original  series  converges  absolutely  and  uniformly  to  some  func- 
tion u(x).  Since  the  n-th  partial  sum  of  the  original  series  is  u  (z),  we  have 

u(x)  =  lim  T\^(x) 

as  n  approaches  infinity. 

All  that  remains  to  "be  shown  is  that  u(x)  satisfies  (II,^).  This  is  easy. 
Consider  the  equation  defining  u^_-(x)  and  let  n  approach  infinity.  We  find  that 

-X  X 

u(x)  =  J     f  (t)[r(x)-r(t)]dt  +  /  q(t)[r(x)-r(t)]q(t)dt, 


which  is  (ll,4).  Consequently  we  have  proved  Theorem  II-l, 

Results  similar  to  those  of  Theorem  II-l  can  he  proved  for  n-th-order  dif- 
ferential equations  or  for  systems  of  first-order  differential  equations,  but  we 
shall  not  discuss  them  here. 
Problem; 

II. 1  Suppose  that  A  is  a  matrix  w:th  piecewise  continuous  elements.  Prove  by 
the  method  of  this  section  that  given  a  constant  vector  U  there  exists  a  vector 
U  such  that  -tt 

dx 
and  such  that  the  value  of  U  for  x  =  0  is  H  ,  For  all  values  of  x  for  which  all 

0 

the  elements  of  A  are  continuous,  the  elements  of  u  are  piecewise  differentiable 
functions  of  x. 

Appendix  III 
Singular  Solutions  of  a  System  of  Equations 

Consider  the  system  of  n  first-order  equations 

g  =  A(x)U. 

where  A  is  a  matrix  whose  elements  are  single-valued  analytic  ftmctions  of  x.  The 
point  X  =  X  is  a  regular  point  of  the  system  if  it  is  a  regular  point  for  all  the 
elements  of  A.  It  is  a  singular  point  of  the  system  if  it  is  a  singularity  of  at 
least  one  of  the  elements  of  A,  It  can  be  shown  that  x  =  x^  is  a  regular  ^ing^^lar 
point  (in  the  same  sense  as  defined  in  the  section  on  singular  points  of  a  differ- 
ential equation)  if  the  singularities  of  the  elements  of  A  are  only  simple  poles. 
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Suppose  for  siarpllcity  that  x  =  0  and  assvune  that  x  is  a  regular  singular 
point  for  the  system.  We  may  expand  each  one  of  the  elements  of  A  in  a  power 
series  around  zero  and  then  collect  terms according  to  powers  of  x.  In  this  way 
we  get  the  following  result : 


|2  =  (A^x'^+A^+A^x  +  A^2+  . . .  )U 


or 


(111.1)  ^  dl  ""  ^W  *  V^  "^  ...)U  =  (A^+  xB)U, 

where  A  ,A  ,A  , ,.,  are  matrices  whose  elements  are  constants  end  B  is  a  matrix 

0   1^ 

whose  elements  are  regular  functions  at  x  =  0. 
Let  us  first  consider  the  simpler  system 

(111. 2)  X  j2  =  A  U  . 
'  dx     0 

If  we  make  the  change  of  independent  variable  x  =  e  ,  (III»2)  becomes 

dt     0 

We  know  that  for  each  eigenvalue  X  'of  the  matrix  A  this  system  has  a  vector 

\t  ° 

solution  with  components  e  ,  If  X  is  an  eigenvalue  of  rank  u,  then  the  system 

also  has  a  solution  with  components  such  as 

Xt  n u  .  _  . u^l 


e 


[t"+  a^t""-^  +  ...  +  aj  . 


Writing  these  results  in  terms  of  x,  we  see  that    (III, 2)  has  a  solution  with 
components  x     or  x  [(log  x)     +  a.  (log  x)        +  ...  +  a^] . 

Let  us  now  return  to  equation   (lll.l).     We  assijae  that   the  vector  U  can 
be  eacpanded  in  a  power  series  as  follows: 

(III.?)  U  =  x°-{X!    +  xU,  +  x\  +   ...). 

ox  ^ 

Here  U  ,  U  ,  U-,...  pre  vectors  whose  components  are  constants.  Note  that 
we  may  assume  U  to  be  not  identically  zero. 

0 

Substituting  the  power  series  for  U  into  (lll.l)  we  get 

ax"v  +    (oEfl)x°^\+   (ccf2)x'^V+   ...  =  x°A  U  +  x°^^(A  U,+  BU  )+x°*^(A  U>BU,  )+< 
0  1  <c  00  olo  o2l 

When  we  equate  the  coefficients  of  corresponding  powers  of  x,  we  obtain  the 
following  set  of  vector  equations  : 
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A  U   =  ctU 
0  0       o 


Vl  =  (a^^V^A 

A^U^  =   (a*2)U2-A^U^-A2U^. 

Each  vector  equation  is  really  a  set  of  n  linear  equations  for  the  n  components 
of  the  corresponding  vector.  Since  the  first  vector  equation  is  homgeneous, 
it  will  have  a  non-zero  solution  if  and  only  if  the  determinant  of  A  -  a,  is  zero. 
This  Implies  that  a  is  an  eigenvalue  of  the  matrix  A  ,  and  then  U^  will  he  an 
eigenvalue  of  A  , 

Let  U  he  any  eigenvalue  of  A  corresponding  to  the  eigenvector  a.  Write 

0  0 

the  second  equation  as  follows: 

(A^-a-l)U^  =  A^U^. 

It  is  clear  that  this  equation  can  he  solved  for  U  if  the  determinant  of  A^-(a  +  l) 
is  not  zero.  Similarly  the  (k  +  l)-th  equation, 

V  "^  -  "-^K  =  "Vk.i  -  Vk-2  -""\\* 

can  he  solved  for  U,  if  the  determinant  of  A  -(a  +  k)  is  not  zero.   In  this  way  we 
may  ohtein  successively  the  vectors  U.,  U-,.,,,  snd  then  the  series  (III. 3)  is 
a  fonnal  solution  of  the  original  system.   It  can  he  proved  that  the  infinite 
series  converges  in  some  neighhorhood  of  the  origin  and  therefore  the  infinite 
series  is  an  actual  solution  of  the  system  in  that  region. 

It  is  clear  that  the  above  argument  is  valid  for  any  eigenvector  of  A^, 
Suppose,  however,  that  A  has  an  eigenvector  of  rank  more  than  one.  What  kind 
of  solution,  if  any,  will  it  correspond  toT  To  make  matters  concrete,  let  us 
assume  that  a  is  an  eigenvalue  of  rank  two  for  A  ,  and  let  us  assume  that 

(Ill.i+)  u  =  z°log  x(U  +  xV^*  x^^   ...)  +  x^(Tq+  xT^+  x^^   ...). 

Substituting  this  expression  in  (lll.l)  and  equating  first  the  coefficients  of 
x"^'^  log  X  and  then  the  coeffients  of  x  '^(n  =  -1,0,1,2,,..),  we  get 
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A  U      =     otU 

0    0  0 


Vl     =      (a+DU^-A^U, 


(III.5) 


A  V       =     aV  +  U 

0    0  0         0 


0* 


Vl     =      (a*lVU^-A^V^ 

AJ2    -     (a  +  2)7^^  U^.  A^V^-  A^V 

•  •  • 

Consider  the  first  equation  of  the  second  set,  namely 

(A^  -  a)V^     =     U  . 

0  0  0 

Since  a  is  an  eigenvalue  of  A  ,  the  determinant  of  A  -  a  vanishes  and  conse^ 

o  0 

quently  we  cannot  conclude  that  this  non-homogeneous  equation  must  have  a  solu- 
tion. However,  we  have  essiisaed  that  a  is  an  eigenvalue  of  rank  two.  let  T 
"be  the  corresponding  eigenvector  of  rank  two;  then  "by  definition 

(A  -  a)^  =  0 


while 


(A^-a)V^  ^  0. 


Denote  the  vector  (A  -  a)V  hy  U  ,  We  see  that  the  first  equation  of  the  first  set 

0      0     0  ^ 

in  (ill, 5)  will  "be  identically  satisfied  since 

(A  -  a)U  =  (A  -  a)\  =  0. 

0      0       0       0 

JTote  that  the  first  set  of  equations  obtained  by  assuming   (III,^)  is  exactly  the 
same  as  the  set  we  obtained  by  assuming   (III,3)«     Let  U(x)  denote  the  solution 
obtained  previously  by  assuming  (III, 3)*   then  from  (III,^)  we  have 

U    =     U(x)log  z  •^  x°'(V  +  xV  +  x^  +  ,,,), 

0         12 

where  V  ,V  "V    ,  satisfies  the  second  set  of  eouations  in  (III,5), 

0   1^ 

Since  we  have  assumed  that  the  determinant  of  A  -(a  +  k)  does  not  vanish,  we 

o 

may  solve  successively  the  remaining  equations  of  the  second  set  to  obtain 
V  (V^,...  .   In  this  way  we  obtain  a  solution  of  the  form  (III,^)  corresponding 
to  any  eigenvector  of  rank  two.  In  the  same  way  can  show  that  if  a  is  an  eigen- 
value of  rank  m  then  -.  t   _  ,  ^-.2 

U  =  U  (log  x)""-^  +  U^  (log  x)"^^  +  ...  -H  U^^ 
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is  a  solution  of  the  system  of  differential  eq-uations.  Here  U  ,  V^,,..,   U  _^ 

are  vectors  whose  components  are  regvUax   functions  of  the  form  x    multiplied 

"by  a  power  series  in  x.     Note  also  that  U  log  x        +  Uj^(log  x)  "  +  .,,  ^^^2 

will  "be  a  solution  of  the  system. 

We  must  now  consider  what  happens  when  "both  a  and  c»fk  are  eigenvalues  of 

A  (k  a  fixed  positive  integer).  In  this  case  the  previous  arguments  break 

o 
down  "because  the  determinant  of  A  -(crt-k)  is  zero.  We  expect  that  the  solution 

0 

will  contain  logarithms  and  we  shall  see  that  this  is  actually  the  case. 
Since  c»fk  is  an  eigenvalue  of  A  ,  there  exists  a  vector 


^o^k  -(k)  ^  ^o^k^^  ^     ^3^,   ^^^j 

0     1      c 


which  is  a  solution  of  the  system. 

Here  W  ,W-,W-,..  are  vectors  whose  components  are  constants.  Note  that  W 
0  1  2  0 

is  an  eigenvector  of  A  corresponding  to  the  eigenvalue  a*k.  We  assume  that 

(111.6)  U  =  x'^CU^H-  x[J^+  ...+  x^"\.i+  3^\+  ...)  +  cx°*^  U^'^^log  X, 

where  c  is  an  arl3itrary  constant.  If  we  substitute  this  expression  in  (III.l) 
we  oh tain  the  equations 

A  U   =  aU 
0  0        o 

Vi  =  («^i)ni-A,u^ 

•  •  • 

Vk  =  (°^^\  -  Vk-1  -  -  -Vo  *  °\- 

The  first  k  of  these  equations  can  he  solved  as  before.     Write  the   (bfl)-th 
eqiiation  as  follows: 

(111.7)  (A^  -  a.k)U^  =  -Vk-r  —  -  Vo  *  °^   • 

Since  the  determinant  of  A  -o-k  vanishes,  this  equation  will  have  a  solution  if 

0 

end  only  if  the  right-hand  side  is  orthogonal  to  all  solutions  of  the  adjoint 

homogeneous  eq-uation. 

Suppose  that  W  is  the  only  eigenvalue  of  A  corresponding  to  the  eigen- 
0  0  » 

value  cc*"k}  then  for  that  same  eigenvalue  there  exists  a  unique  eigenvector  W 

of  the  adjoint  matrix  A  .  Now,  if  we  choose  c  so  that 

c(W*,  W^)  =  (\.h\^l^   "'*W 
then  the  right-hsjidside  of  (ill. 7)  will  he  orthogonal  to  the  solutions  of  the 
adjoint  homogeneous  equation  and  therefore  there  will  exist  a  vector  Uj^  satisfy- 
ing (III. 7). 


-  I8i^  - 

The  remaining  equations  for  U   ,  Uj_2»  •••  '^an  be   solved  as  "before,  since 
■by  assTomption  det[A  -dr-lz-mj    (m  =  1,2,,,,)  does  not  vanish.  In  this  way  we  have 
obtained  a  series  of  the  form  (ill, 6)  which  is  a  formal  solution  of  the  differ- 
ential system.  Again,  it  can  be  proved  that  this  series  will  converge  in  some 
neighborhood  of  the  origin  and  therefore  it  will  act\ially  be  a  solution. 

We  combine  our  resiolts  in  the  following  theorem: 
Theorem I I I-l,  Consider  the  following  system  of  first-order  differential  ecuations; 

(III.8)  xg  =  (A^  +  iB)U. 

where  A  is  a  matrix  whose  elements  are  constants  and  B  is  a  matrix  whose  elements 
0 ' — — ■ ' — — — — 

are  re/ralar  functions  of  x  in  the  neighborhood  of  x  =  0,  Let  a  be  an  eigenvalue 

of  A  such  that  a  -<•  m.  (m  any  Dositive  Integer)  is  not  an  eigenvalue  of  A  ;  then. 
0  0 

if  U  is  any  eigenvector  of  A  corresponding  to  the  eigenvalue  a.  (I II, 8)  has  a, 
solution  n(x)t 

(III.8)  U(x)  =  x°-{-U+   xU.  +  ...). 

o    1 

If  7  is  an  eigenvector  of  rank  p  correspondine:  to  the  eigenvalue  a  (where  a 

is  such  that  a-  •'•m  is  not  an  eigenvalue),  then  there  exists  a  solution  of  (III. 3) 

of  the  form 

U^(x)(log  x)°-^  +  U^(x)(log  x)°^2  ^  ^^^  ^  U^^(x)(x), 

where 

U^(x)(log  x)"^2  ■••  U^(x)(log  x)"^^  +  ,.,  +  U^2^^^^^^ 

is  also  a  solution  of  (III, 8).  and  where 

•  ^„^t(^)  =  3C°(7  +  xV  +  x^_  +  .,.). 
m— 1  0    1     c 

Let  g  be  an  eigenvalue  of  A  such  that  a  +  k  (k  a  positive  integer)  is  also  a 
eigenvalue.  If  U  Is  a  eigenvector  of  A  correspondige:  to  the  eigenvalue  a,  then 
there  exists  a  solution  of  (III. 3)  of  the  form 

x*^(U  +  xU,  +  ,,. )  +  x°^^  U,  log  X, 

0     1  ^ 

where  x    U  is  also  a  solution  of  ,  (111,8). 
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